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Abstract 

We investigate the phase-space of a flat FRW universe including both a scalar field, 
4>, coupled to matter, and radiation. The model is inspired in scalar-tensor theories 
of gravity, and thus, related with F(R) theories through conformal transformation. 
The aim of the chapter is to extent several results to the more realistic situation when 
radiation is included in the cosmic budget particularly for studying the early time dy- 
namics. Under mild conditions on the potential we prove that the equilibrium points 
corresponding to the non-negative local minima for V(<j>) are asymptotically stable. 
Normal forms are employed to obtain approximated solutions associated to the in- 
flection points and the strict degenerate local minimum of the potential. We prove 
for arbitrary potentials and arbitrary coupling functions x(4>)i °f appropriate differen- 
tiable class, that the scalar field almost always diverges into the past. It is designed 
a dynamical system adequate to studying the stability of the critical points in the 
limit \cf>\ — > oo. We obtain there: radiation-dominated cosmological solutions; power- 
law scalar-field dominated inflationary cosmological solutions; matter-kinetic-potential 
scaling solutions and radiation-kinetic-potential scaling solutions. Using the mathemat- 
ical apparatus developed here, we investigate the important examples of higher order 
gravity theories F(R) = R + aR 2 (quadratic gravity) and F (R) = R n . We illustrated 
both analytically and numerically our principal results. In the case of quadratic gravity 
we prove, by an explicit computation of the center manifold, that the equilibrium point 
corresponding to de Sitter solution is locally asymptotically unstable (saddle point). 
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1 Introduction 

Modern Cosmology is a broad and promising area of research in applied mathematics and 
applied physics based on the observed data available in the modern astronomical literature. 
One of the most celebrated discoveries of observational cosmology is that the observable 
universe is now known to be accelerating [31 131 IH [S] , and this feature led physicists to follow 
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two directions in order to explain it. The first is to introduce the concept of dark energy 
(see the reviews [51 [71 [HI HI E3 an d references therein) in the right-hand-side of the field 
equations, which could either be the simple cosmological constant A or, one or several 
scalar fields [HI [TBI H7J EH E3 ■ The second is looking for alternative models [20 1 12T 1 122 ] . 
Alternative approaches to dark energy are the so-called Extended Theory of Gravitation 
(ETG) and, in particular, higher-order theories of gravity (HOG) G3j Ell [2j EBJ [27J EH 

[^02312113211121132113312^ Such an 

approach can still be in the spirit of General Relativity Theory (GRT) since the only request 
the Hilbert-Einstein action should be generalized (by including non-linear terms in the 
Ricci curvature R and/or involving combinations of derivatives of R (48j [49j [50l [5T] ) asking 
for a gravitational interaction acting, in principle, in different ways in both cosmological 
[25] 126} 152"! l53l 154} 1551 [56] and astrophysical [52j [57] scales. In this case the field equations 
can be recast in a way that the higher order corrections are written as an energy-momentum 
tensor of geometrical origin describing an "effective" source term on the right hand side of 
the standard Einstein field equations [25j [26] . These models have been studied from the 
dynamical systems viewpoint in gjj OH ESI M ED ES El E3 El E3 E3 ED El EH] ■ On the 
other hand, scalar fields have played an essential role as models for the early-time universe. 
In the inflationary universe scenarios (mainly based on GRT) matter is modeled, usually, as 
a scalar field, <fi, with potential V(<fi), which must meet the requirements necessary to lead 
to the early-time accelerating expansion [TCJ [ZD (72J [731 171] • If the potential is constant, 
i.e., if V{4>) = Vq, space-time is de Sitter and expansion is exponential. If the potential 
is exponential, i.e., V(4>) = Vq exp[— \4>], we get an inflationary powerlaw solution [751 176] , 
"Extended" inflation models [77J E3 HSl 1HD] , on the other hand, use the Brans-Dicke theory 
(BDT) [81] as the correct theory of gravity, and in this case the vacuum energy leads directly 
to a powerlaw solution |82j while the exponential expansion can be obtained if a cosmological 
constant is explicitly inserted into the field equations [77j [83j |84] . 

Several gravity theories consider multiple scalar fields with exponential potential, par- 
ticularly assisted inflation scenarios [85l [86l HTJ [88l [89], quintom dark energy paradigm 
[90l 191] l92l 193] and others. Also, have been considered positive and negative exponen- 
tial potentials [M], single exponential and double exponential [HI (SBJ H7J [§B] , etc. Other 
generalizations with multiple scalar fields are available [99l II 00] . 

In BDT, a scalar field, \i ac ts as the source for the gravitational coupling with a varying 
Newtonian 'constant' G ~ X 1 - BDT is the first propotype of the so-called Scalar- Tensor 
theories (STT) of gravity jTJJD EH EH EH EH EE]. It is worthy to mention that 
BDT survive several observational tests including Solar System tests |107j and Big-Bang 
nucleosynthesis constraints 110811109] , More general STT with a non-constant BD parameter 
w(x), and non-zero self-interaction potential V(x)t have been formulated, and also survive 
astrophysical tests [TTOl ITT p lTT2 ] . 

To our knowledge, the dynamical behavior of space-times based on GRT is so far known 
for a large variety of models with scalar fields with non-negative potential |113[ I114[ I115[ 
11161 11171 1118) . In reference [118] , have been extended many of the results obtained in |114j 
considering arbitrary potentials. In [113 it has been shown that for a large class of FRW 
cosmologies with scalar fields with arbitrary potential, the past attractor is a family of 
solutions in one-to-one correspondence with exactly integrable cosmologies with a massless 
scalar field. This result has been extended somewhat in [1] to FRW cosmologies based 
on STTs. In this reference was investigated a general model of coupled dark energy with 
arbitrary potential V and coupling function \- It was proved there, by using dynamical 
systems techniques that if the potential and the coupling function are sufficiently smooth 
functions; the scalar field almost always diverges into the past. Under some regularity 

1 This choice is seriously plagued by the well known coincidence and fine tuning problems [111 1121 1131 1 141 . 
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conditions for the potential and for the coupling function in that limit, it was constructed 
a dynamical system well suited to investigate the dynamics near the initial singularity. 
The critical points therein were investigated and the cosmological solutions associated to 
them were characterized. There was presented asymptotic expansions for the cosmological 
solutions near the initial space-time singularity, which extend previous results of |113) . On 
the other hand, in |119j it was investigated flat and negatively curved Friedmann-Robertson- 
Walkcr (FRW) models with a perfect fluid matter source and a scalar field arising in the 
conformal frame of F(R) theories nonminimally coupled to matter. It was proved there 
that, for a general class of potentials V, the equilibrium corresponding to non- negative local 
minima for V are asymptotically stable, as well as horizontal asymptotes approached from 
above by V. For a nondegenerate minimum of the potential with zero critical value they 
prove in detail that if 7 > 1, then there is a transfer of energy from the fluid to the scalar 
field and the later eventually dominates in a generic way. As we will see in next sections 
the results in 119J and in [1 can be obtained by investigating a general class of models 
containing both STTs and F(R) gravity. 

The aim of the chapter is to extent several results in |113llll4l [Tl lll9j to the more realistic 
situation when radiation is included in the cosmic budget (particularly for investigating the 
early-time dynamics). We will focus mainly in a particular era of the universe where matter 
and radiation coexisted. Otherwise the inclusion of radiation complicates the study in an 
unnecessary manner, since, assuming a perfect barotropic fluid with an arbitrary barotropic 
index 7, for 7 = 4/3, this matter source corresponds to radiation. Thus we will consider both 
ordinary matter described by a perfect fluid with equation of state p = (7— l)p (coupled to a 
scalar field) and radiation with energy density p r . We are interested in investigate all possible 
scaling solutions in this regime. Although we are mainly interest in describing the early 
time dynamics of our model, for completeness we will focus also in the late-time dynamics. 
As in [119] we obtain for flat FRW models sufficient conditions under the potential, to 
establish the asymptotic stability of the non- negative local minima for V{<p). Normal forms 
are employed to obtain approximated solutions associated to the local degenerated minimum 
and the inflection points of the potential. We prove for arbitrary potentials and arbitrary 
coupling functions x(0)i °f appropriate differentiable class, that the scalar field almost 
always diverges into the past generalizing the results in [1131 fTj). It is designed a dynamical 
system adequate to studying the stability of the critical points in the limit \<j>\ — > 00. We 
obtain there: radiation-dominated cosmological solutions; power-law scalar-field dominated 
inflationary cosmological solutions; matter-kinetic-potential scaling solutions and radiation- 
kinctic-potential scaling solutions. It is discussed, by means of several worked examples, 
the link between our results and the results obtained for specific F(R) frameworks by using 
appropriated conformal transformations. We illustrated both analytically and numerically 
our principal results. Particularly, we investigate the important examples of higher order 
gravity theories F(R) = R + aR 2 (quadratic gravity) and F(R) = R n . In the case of 
quadratic gravity we prove, by an explicit computation of the center manifold, that the 
equilibrium point corresponding to de Sitter solution is locally asymptotically unstable 
(saddle point). This result complements the result of the proposition discussed in [115] p. 
5, where it was proved the local asymptotic instability of the de Sitter universe for positively 
curved FRW models with a perfect fluid matter source and a scalar field which arises in the 
conformal frame of the R + aR 2 theory. Finally, we investigate a general class of potentials 
containing the cases investigated in [165, 166 . In order to provide a numerical elaboration 
for our analytical results for this class of models, we re-examine the model with power-law 
coupling and Albrecht-Skordis potential V{4>) — e~^{A + (tj>— B) 2 ) investigated in [1] in 
presence of radiation. 
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2 Brief review on dynamical systems analysis of cosmo- 
logical models 

Qualitative methods have been proved to be a powerful scheme for investigating the physical 
behavior of cosmological models. It has been used three different approaches: approximation 
by parts, Hamiltonian methods, and dynamical systems methods [120 . In the third case 
the Einstein's field equations of Bianchi's cosmologies and its isotropic subclass (FLRW 
models), can be written as an autonomous system of first-order differential equations whose 
solution curves partitioned to K™ in orbits, defining a dynamical system in R™. In the 
general case, the elements of the phase space partition (i.e., critical points, invariant sets, 
etc.) can be listed and described. This study consists of several steps: determining critical 
points, the linearization in a neighborhood of them, the search for the eigenvalues of the 
associated Jacobian matrix, checking the stability conditions in a neighborhood of the critical 
points, the finding of the stability and instability sets and the determination of the basin of 
attraction, etc. In some occasions, in order to do that, it is needed to simplify a dynamical 
system. Two approaches are applied to this objective: one, reduce the dimensionality of 
the system and two, eliminate the nonlinearity. Two rigorous mathematical techniques that 
allow substantial progress along both lines are center manifold theory and the method of 
normal forms. We will use both approaches in our analysis. We submit the reader to sections 
I2.1.5l and l2.1.4l for a summary about such techniques. 

The most general result to be applied in order to determine the asymptotic stability of a 
equilibrium point, a, is Liapunov's stability theorem. Liapunov's stability method provides 
information not only about the asymptotic stability of a given equilibrium point but also 
about its basin of attraction. This cannot be obtained by the usual methods found in the 
literature, such as linear stability analysis or first-order perturbation techniques. Moreover, 
Liapunov's method is also applicable to non-autonomous systems. To our knowledge, there 
are few works that use Liapunov's method in cosmology [1211 11221 11231 1 124=1 1125) . In |124j 
it is investigated the general asymptotic behavior of Friedman-Robertson- Walker (FRW) 
models with an inflaton field, scalar-tensor FRW cosmological models and diagonal Bianchi- 
IX models by means of Liapunov's method. In [125j it is investigated the stability of 
isotropic cosmological solutions for two-field models in the Bianchi I metric. The author 
proved that the conditions sufficient for the Lyapunov stability in the FRW metric provide 
the stability with respect to anisotropic perturbations in the Bianchi I metric and with 
respect to the cold dark matter energy density fluctuations. Sufficient conditions for the 
Liapunov's stability of the isotropic fixed points of the system of the Einstein equations are 
also found (these conditions coincided with the previously obtained in [93] for the quintom 
paradigm without using Liapunov's technique). To apply Liapunov's stability method it 
is required the construction of the so-called strict Liapunov's function, i.e., a C 1 function 
V : U C R n ->• K defined in a neighborhood U of a such that V(a) = 0, V(x) >0,i/a 
and V(x) < (< 0) in U \ {a}. The construction of such V is laborious, and sometimes 
impossible. One alternative way is to follow point (5) in section l2~2l 

For the investigation of hyperbolic equilibrium points of autonomous vector fields we 
can use Hartman-Grobman's theorem (theorem 19.12.6 en [126] p. 350) which allows for 
analyzing the stability of an equilibrium point from the linearized system around it. For 
isolated nonhyperbolic equilibrium points we can use normal forms theorem (theorem 2.3.1 
in [127j ). which contains Hartman-Grobman's theorem as particular case. The aim of the 
normal form calculation is to construct a sequence of non-linear transformations which 
successively remove the non-resonant terms in the vector field of order r in the Taylor's 
expansion, X r , starting from r = 2. Normal forms (NF) theory has been used in the context 
of cosmological models in order to get useful information about isolated nonhyperbolic 
critical points. In [128] was investigated Hofava-Lifshitz cosmology from the dynamical 
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system view point. There was proved the stability of dc Sitter solution (corresponding to 
a non-hyperbolic critical point) in the case where the detailed-balance condition is relaxed, 
using NF expansions. In |129j were investigated closed isotropic models in second order 
gravity. There the normal form of the dynamical system has periodic solutions for a large 
set of initial conditions. This implies that an initially expanding closed isotropic universe 
may exhibit oscillatory behavior. 

The more relevant concepts of qualitative theory are the concept of flow, and the concept 
of invariant manifold. The invariant manifold theorem (theorem 3.2.1 en |126) h that claims 
for the existence of local stable and unstable manifolds (under suitable conditions for the 
vector field), only allows to obtain partial information about the stability of equilibrium 
points and does not gives a method for determining the stable and unstable manifolds. 
Sometimes it is required to consider higher order terms in the Taylor's expansion of the 
vector field (e.g. normal forms). For the investigation of the asymptotic states of the system 
the appropriated concepts are a and w-limit sets of a; £ M. n (definition 8.1.2 en [126] p. 105). 
To characterize these invariant sets one can use the LaSalle's Invariance Principle or the 
Monotonicity Principle f |120j . p. 103; |130j p. 536). To apply the Monotonicity Principle 
it is required the construction of a monotonic function. In some cases a monotonic function 
is suggested by the form of the differential equation (see equation (fTF]) below) and in some 
cases by the Hamiltonian formulation of the field equations. In |131) is given a prescription 
of how to find monotonic functions for Bianchi cosmologies using Hamiltonian techniques, 
merely as a convenient tool for an intermediate step; the final results are described in 
terms of scale-automorphism invariant Hubble-normalized reduced state vector, which is 
independent of a Hamiltonian formulation. Also, one can use the Poincare-Bendixson ( [132] 
p. 22, theorem 2 |133j p. 6, [134] ) theorem and its collorary to distinguishing among all 
possible w-limit sets in the plane. From the Poincare-Bcndixon Theorem follows that any 
compact asymptotic set is one of the following; 1.) a critical point, 2.) a periodic orbit, 
3.) the union of critical points and heteroclinic or homoclinic orbits. As a consequence if 
the existence of a closed (i.e., periodic, heteroclinic or homoclinic) orbit can be ruled out 
it follows that all asymptotic behavior is located at a critical point. To ruled out a closed 
orbits for two-dimensional systems we can use the Dulac's criterion (theorem 3 [133 p. 6, 
se also [120] . p. 94, and |132j ). It requires the construction of a Dulac's function. A Dulac's 
function, B, is a C 1 function defined in a simply connected open subset D C R 2 such that 
V W) = kWi) + £^( B M > 0, or (< 0) for all x G D (see [133]). 

Finally, in order to the chapter to be self-contained, let us formulate (without proof) 
some of the results outlined above that we shall use in our proofs. 

2.1 Some terminology and results from the Dynamical Systems 
Theory 

In the following 0(r, x) denotes the flow generated by the vector field (or differential equa- 
tion) 

x'(r) = f(x(r)), x(r) G R", (1) 
where the prime denote derivative with respect to t. 

2.1.1 Limit sets 

Definition 1 (definition 8.1.1, [126J p. 104) A point xo G K n is called an u-limit point 
of x G R™, denoted w(x), if there exists a sequence {t^}, Tj — > oo such that (f>(ri,x) — > xq. 
a-limits are defined similarly by taking a sequence {t^}, Tj — > — oo. 
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Definition 2 (definition 8.1.2, [126J p. 105) The set of all oj- limit points of a flow or 
map is called a u-limit set. The a-limit is similarly defined. 

Proposition 1 (proposition 8.1.3 [126J, p. 105) Let </>(r, •) be a flow generated by a 
vector field and let M be a positively invariant compact set for this flow (see definition 3.0.3 
p. 28 JIWI ). Then for p € M, we have 

i) w(p) + 
ii) w(p) is closed 

Hi) w(p) is invariant under the flow, i.e., w(p) is a union of orbits, 
iv) U)(p) is connected. 

2.1.2 Monotone functions and Monotonicity Principle 

Definition 3 (definition 4.8 [120J, p. 93) Let <j)(r, x) be a flow on R", let S be an in- 
variant set of 4>(t, x) and let Z : S — > R be a continuous function. Z is monotonic decreasing 
(increasing) function for the flow <P(t, x) means that for all x € S, Z((f>(r, x)) is a monotonic 
decreasing (increasing) function of t. 

Proposition 2 (Proposition 4.1, [120J, p. 92) Consider a differential equation (JXJ) with 
flow 0(r,x). Let Z : R" —> R be a C 1 (R™) function which satisfies Z' = aZ, where 
a : R ra — >■ R is a continuous function. Then, the subsets o/M" defined by Z > 0, Z — 0, or 
Z < are invariant sets for </)(r, x). 

Theorem 1 (Monotonicity Principle) Let 4>(t, x) be a flow on R n with S an invariant 
set. Let Z : S — > R be a C 1 (R") function whose range is the interval (a, b) where a G 
R U {— oo}, b G R U {+oo}, and a < b. If Z is decreasing on orbits in S, then for all x € S, 
w(x) c{seS \ S\lim y ^ s Z(y) ^ b} and a(x) C {s E S \ S\lim y ^ s Z(y) ^ a}. 

2.1.3 Results for planar systems 

Theorem 2 (Poincare-Bendixon Theorem) Consider the differential equation (p} on 
R 2 , with f G C 2 , and suppose that there are at most a finite number of singular points (i.e., 
no non-isolated singular points). Then any compact asymptotic set is one of the following; 
1. a singular point, 2. a periodic orbit, 3. the union of singular points and heteroclinic 
or homoclinic orbits. 

Theorem 3 (Dulac's Criterion) If D C R 2 is a simply connected open set and B is a 
Dulac's function on D, then the differential equation ([T]) on I 2 , with f G C 1 has no periodic 
(or closed) orbit which is contained in D. 

2.1.4 Normal Forms 

In this section we offer the main techniques for the construction of normal forms for vector 
fields in R" . We follow the approach in [127j . 

Let X : R™ — > R" be a smooth vector field satisfying X(0) = 0. We can formally 
construct the Taylor expansion of x around 0, namely, X = Xi +X2 + . . . + Xfc + 0(|x| fc+1 ), 
where X r G H r , the real vector space of vector fields whose components are homogeneous 
polynomials of degree r. For r — 1 to k we write 
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X r (x) — X)mi = l ■ • ■ Em„=l Ej=l ■^■m,jX m ej, 

m = r, (2) 
Observe that Xi = DX(0)x = Ax, i.e., the Jacobian matrix. 

The aim of the normal form calculation is to construct a sequence of transformations 
which successively remove the non-linear term X r , starting from r = 2. 
The transformation themselves are of the form 

x = y + h r (y), (3) 

where h r £ H r , r > 2. 

The effect of ((3|) in Xi is as follows [ 1 2 7 j : Observe that x = 0(|y|). Then, the inverse of 
([3]) takes the form 

y = x-h r (x)+0(|xr +1 ). (4) 
By applying total derivatives in both sides, and assuming x' = Ax + X r (x), we find 

y' = Ay - L A h r (y) + X r (y) + 0(|y|' r+1 ) (5) 

where La is the linear operator that assigns to h(y) £ H r the Lie bracket of the vector 
fields Ay and h(y): 

L A : H r -> H r 

h ^L A h(y) =Dh(y)Ay-Ah(y). (6) 

Both La and X r £ H r , so that the deviation of the right-hand side of (|5|) from Ay has 
no terms of order less than r in |y|. This means that if X is such that X2 = . . . X r _i = 0, 
they will remain zero under the transformation ([3]). This makes clear how we may be able 
to remove X r from a suitable choice of h r . 

The proposition 2.3.2 in |127j states that if the inverse of L A exists, the differential 
equation 

x' = Ax + X r (x)+0(|x| r+1 ) (7) 
with X r £ H r ', can be transformed to 

y' = Ay + 0(|yr +1 ) (8) 

by the transformation Q where 

h r (y) =L A - 1 X,.(y) (9) 

The equation 

L A h r (y) = X r (y) (10) 

is named the homological equation. 

If A has distinct eigenvalues Aj, i — 1, 2, 3, its eigenvectors form a basis of K". Relative to 
this eigenbasis, A is diagonal. It can be proved (see proof in |127j ) that La has eigenvalues 
A m j = m • A — Xi — ^2 ■ rrijXj — \i with associated eigenvectors x m ei. The operator, La ~ , 
exists if and only if the A m i ^ 0, for every allowed m and i = 1 . . . r. 

If we were able to remove all the nonlinear terms in this way, then the vector field can 
be reduced to its linear part 

x' = X(x) -> y' = Ay. 
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Unfortunately, not all the higher order terms vanishes by applying these transformations. 
It is the case if resonance occurs. 

The n-tuple of eigenvalues A = (Ai, . . . , A„) T is resonant of order r (sec definition 2.3.1 
in |127| ) if there exist some m = (mi, m^, ■ ■ ■ rn n ) T (a n-tuple of non-negative integers) with 
mi + ni2 + . . . m n = r and some i = 1 . . . n such that Ai = m • A, i.e., if A m i = for some 
m and some i. 

If there is no resonant eigenvalues, and provided they are different, we can use the 
eigenvectors of A as a basis for H r . Then, we can write h r as 

h r (x) = ^2 /i m ,iX m ej 

m,i,Y2 m j — r 

and any vector field X E H r as 

X(x) = J2 X m.X m ei 

m,i, ^2 rnj=r 

where m = (mi,i7i2, . . . , m n ) T , x m = a;™ 1 x™ 2 . . . x™ n and e*, i = 1, . . . n stands for the 
canonical basis in R™ . If the eigenvalues of A are not resonant of order r, then 

hm,i — ^m,z/Am,2- 

This gives h r explicitly in terms of X r . 

In case that resonance occurs, we proceed as follows. If A can diagonalized, then the 
eigenvectors of La form a basis of H r . The subset of eigenvectors of La with non-zero eigen- 
values then form a basis of the image, B r , of H r under La- It follows that the component 
of X r in B r can be expanded in terms of these eigenvectors and h r chosen such that 

hm,i — ^m,i/Am,z- 

to ensure the removal of these terms. The component, w r , of X r lying in the complementary 
subspace, G r , of B r in H r will be unchanged by the transformations x = y + h. r (y) obtained 
from B r . 
Since 

X r (y + h r+fe (y)) = X r (y) + 0(| y r +fc+1 ), r > 2, k = 1, 2, . . . , 

these terms are not changed by subsequent transformations to remove non-resonant terms 
of higher order. 

The above facts are expressed in 

Theorem 4 (theorem 2.3.1 in |127j ) Given a smooth vector field X(x) onR™ withX(0) = 0, 
there is a polynomial transformation to new coordinates, y, such that the differential equa- 
tion x' = X(x) takes the form y' = Jy+X^i w r(y) + 0(\y\ N+1 ), where J is the real Jordan 
form of A = DX(0) and w r 6 G r , a complementary subspace of H r on B r = \^^(H r ). 

2.1.5 Center Manifold theory 

In this section we offer the main techniques for the construction of center manifolds for 
vector fields in K™. We follow the approach in [126j chapter 18. 
The setup is as follows. We consider vector fields in the form 

x' = Ax + f(x,y), 

y' = Bx + g(x,y), (x,y)GM c xR s , (11) 
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where 



f(0,0) = 0,Df(0,0) = 0, 
g(0,0) =0,Dg(0,0) = 0. 



(12) 



In the above, A is a c x c matrix having eigenvalues with zero real parts, B is an s x s 
matrix having eigenvalues with negative real parts, and f and g are C r functions (r > 2). 

Definition 4 (Center Manifold) An invariant manifold will be called a center manifold 
for ([lip if it can locally be represented as follows 

W c (0) = {(x, y) e K c x M s : y = h (x) , |x| < 5} ; h (0) = 0, Dh (0) = 0, 

for 6 sufficiently small (cf. JJgtf p. 246, fWSf.v. 155). 

The conditions h (0) = 0, £>h(0) = imply that W c (0) is tangent to E c at (x,y) = 
(0,0), where E c is the generalized eigenspace whose corresponding eigenvalues have zero 
real parts. The following three theorems (see theorems 18.1.2, 18.1.3 and 18.1.4 in |126] 
p. 245-248) are the main results to the treatment of center manifolds. The first two are 
existence and stability theorems of the center manifold for (|11[) at the origin. The third 
theorem allows to compute the center manifold to any desired degree accuracy by using 
Taylor series to solve a quasilinear partial differential equation that h (x) must satisfy. The 
proof of those results is given in [135] . 

Theorem 5 (Existence) There exists a C r center manifold for (|11[) . The dynamics of 
(JTTJ) restricted to the center manifold is, for u sufficiently small, given by the following 
c-dimensional vector field 



The next results implies that the dynamics of (|13|) near u = determine the dynamics 
of (HI]) near (x,y) = (0,0) (see also Theorem 3.2.2 in [164] ). 

Theorem 6 (Stability) i) Suppose the zero solution of (jl3[) is stable (asymptotically sta- 
ble) (unstable); then the zero solution of (|11|) is also stable (asymptotically stable) (unstable) . 
Then if (x(t), y(r)) is a solution of (jlip with (x(0),y(0)) sufficiently small, then there is 
a solution u(t) of (|13|) such that, as r — » oo 



where r > is a constant. 

Dynamics Captured by the center manifold 

Stated in words, this theorem says that for initial conditions of the full system sufficiently 
close to the origin, trajectories through them asymptotically approach a trajectory on the 
center manifold. In particular, equilibrium points sufficiently close to the origin, sufficiently 
small amplitude periodic orbits, as well as small homoclinic and heteroclinic orbits are 
contained in the center manifold. 

The obvious question now is how to compute the center manifold so that we can use 
the result of theorem |6p To answer this question we will derive an equation that h(x) must 
satisfy in order to its graph to be a center manifold for (|11[) . 



u' = Au + f (u, h (u)) , u£l c . 



(13) 



x(r)=u(r)+0( e - rt ), 
x(r)=h(u(r))+0( e - rt ), 
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Suppose we have a center manifold 

W c (0) = {(x, y) G K c x R s : y = h (x) , |x| < 6} ; h (0) = 0, Dh (0) = 0, 

with S sufficiently small. Using the invariance of W c (0) under the dynamics of fjlljt . we 
derive a quasilinear partial differential equation that h (x) must satisfy. This is done as 
follows: 

1. The (x, y) coordinates of any point on W° (0) must satisfy 

y = h(x) (14) 

2. Differentiating (TT41 with respect to time implies that the (x',y') coordinates of any 
point on W c (0) must satisfy 

y' = Dh(x)x' (15) 

3. Any point in W c (0) obey the dynamics generated by (fTTj) . Therefore substuting 

x' = Ax + f (x,h(x)) , 
y' = Bh(x)+g(x,h(x)) 

into (|T5|) gives 

N (h(x)) = Dh(x) [Ax + f (x, h(x))] - Bh(x) - g (x, h(x)) = 0. (16) 

Equation (fT!)|) is a quasilinear partial differential that h(x) must satisfy in order for its 
graph to be an invariant center manifold. To find the center manifold, all we need to do is 
solve HU). 

Unfortunately, it is probably more difficult to solve (|16[) than our original problem; 
however the following theorem give us a method for computing an approximated solution 
of (fl6|) to any desired degree of accuracy. 

Theorem 7 (Approximation) Let <& : R c — > R s be a C 1 mapping with 4*(0) = and 
D*(0) = such that M (*(x)) = 0{\x\ q ) asx^O for some q > 1. Then, |h(x) -*(x)| = 
0(\x\ q ) asx^O 

This theorem allows us to compute the center manifold to any desired degree of accuracy 
by solving (JTHJ) to the same degree of accuracy. For this task power series expansions will 
work nicely. Let us consider a concrete example further in section [4.2.1 1 



2.2 Standard procedure to analyze the properties of the flow of a 
vector field 

Given a cosmological dynamical system determined by the differential equation 



^ = f(y),yeR n , (17) 

g(y) = o, (is) 

the standard procedure to analyze the properties of the flow generated by (|17|) subject to 
the constraint(s) (fT8|) (see, for example, the reference [120 ) is the following: 



1. Determine whether the state space, as defined by (|18|) . is compact 
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2. Identify the lower-dimensional invariant sets, which contains the orbits of more special 
classes of models with additional symmetries. 

3. Find all the equilibrium points and analyze their local stability. Where possible iden- 
tify the stable and unstable manifolds of the equilibrium points, which may coincide 
with some of the invariant sets in point (2). 

4. Find Dulac's functions or monotone functions in various invariant sets where possible. 

5. Investigate any bifurcation that occur as the equation of state parameter 7 (or any 
other parameters) varies. The bifurcations are associated with changes in the local 
stability of the equilibrium points. 

6. Having all the information in the points (l)-(5) one can hope to formulate precise 
conjectures about the asymptotic evolution, by identifying the past and the future 
attractors. The past attractor will describe the evolution of a typical universe near 
the initial singularity while the future attractor will play the same role at late times. 
The monotone functions in point (4) above, in conjunction with theorems of dynamical 
systems theory, may enable some of the conjectures to be proved. 

7. Knowing the stable and unstable manifolds of the equilibrium points it is possible 
to construct all possible heteroclinic sequences that join the past attractor, thereby 
gaining insight into the intermediate evolution of cosmological models. 

3 Flat FRW models with both a scalar field coupled to 
matter and radiation 

There exist a number of physical theories which predict the presence of a scalar field coupled 
to matter. For example, in string theory the dilaton field is generally coupled to matter 
|136j . Nonminimally coupling occurs also in STT of gravity |137( 1138) , in HOG theories [27] 
and in models of chameleon gravity [139] . Coupled quintessence was investigated also in 
|140i 11411 1142] by using dynamical systems techniques. The cosmological dynamics of scalar- 
tensor gravity have been investigated in [5911143] . Phenomenological coupling functions were 
studied for instance in [144] which can describe either the decay of dark matter into radiation, 
the decay of the curvaton field into radiation or the decay of dark matter into dark energy 
(see section III of [144] for more information and for useful references). In the reference [143] . 
the authors construct a family of viable scalar-tensor models of dark energy (which includes 
pure F(R) theories and quintessence). By investigating a phase space the authors obtain 
that the model posses a phase of late-time acceleration preceded by a standard matter era, 
while at the same time satisfying the local gravity constraints (LGC). 

The aim of the section is to extent several results in [11311114] [T| I119] to the more realistic 
situation when radiation is included in the cosmic budget (particularly for the study of the 
early time dynamics). We will focus mainly in a particular era of the universe where matter 
and radiation coexisted. We are interested in investigate all possible scaling solutions in 
this regime. 

3.1 The action 

The action for a general class of STT, written in the so-called Einstein frame (EF), is given 
by [145]: 

S E F= J m ^x/Mj^-^V^^y^+xW^^atter^V^X^-^a/i)! (19) 
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In this equation R is the curvature scalar, (f> is the a scalar field, related via conformal 
transformations with the dilaton field, x-0 V(4>) is the quintessence self-interaction potential, 
X(4>)~ 2 is the coupling function, £ ma tter is the matter Lagrangian, p. is a collective name 
for the matter degrees of freedom. The energy-momentum tensor of background matter is 
defined by 

2 5 ( 

T aP = - ^ {^X-'Am.Vm,*- 1 ^)} • (20) 

In the STT written in the EF as (TT9"1) . the conservation equations for matter/energy are 
expressed in 

Qp = V T Q(3 = --T V^, 

where T = T" denotes the trace of the energy-momentum tensor of background matter. For 
HOG theories derived from Lagrangians of the form 

L = -F (R) + Cmatter(fJ; V/i, g a p), (21) 

it is well known that under the conformal transformation, g^ = F' (i?) g^ , the field equa- 
tions reduce to the Einstein field equations with a scalar field 4> as an additional matter 
source, where 

<t>=sJ\\nF'{R). (22) 

Assuming that (|2"2")l can be solved for R to obtain a function R ((f)) , the potential of the 
scalar field is given by 

V(R(<f>)) = —^s(RF'-f), (23) 
and quadratic gravity, F(R) = R + aR 2 with the potential 

is a typical example. The restrictions on the potential in the papers [116, 1561 H57[ 1117] 
were used in |158j to impose conditions on the function / (i?) with corresponding potential 
(1231) . The conformal equivalence can be formally obtained by conformally transforming the 
Lagrangian (|2Tj) and the resulting action becomes 159 , 

S = J d 4 xy/^ | Ir - (V0) 2 -V{4>) + e- 2 ^C matt er (fh Vp, e" I 

It is easy to note that the model arising from the action (pM]l can be obtained from (H 
with the choice x{<t>) = ev^*. 



(24) 



3.2 The field equations 

In this section we consider a phenomenological model inspired in the action (TT91 for FRW 
space-times with flat spatial slices, modelled by the metric: 

ds 2 = -dt 2 + a(i) 2 (dr 2 + r 2 (d9 2 + sin 2 9dip 2 )) . (25) 

2 For a discussion about the regularity of the conformal transformation, or the equivalence issue of the 
two frames, see for example [1461 [1471 ITl8l [1491 liliOl 11511 [1521 [T53l [T54l [T55] and references therein. 
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We use a system of units in which 8irG = c = h = 1. We assume that the energy- momentum 
tensor ([2l7| is in the form of a perfect fluid 



Tp = diag (-p,p,p,p) , 

where p and p are respectively the isotropic energy density and the isotropic pressure (con- 
sistently with FRW metric, pressure is necessarily isotropic |160j ). For simplicity we will 
assume a barotropic equation of state p = (7 — l)p. We include radiation in the cosmic 
budget since it is an important matter source in the early universe. 

The cosmological equations for flat FRW models with both a scalar field coupled to 
matter and radiation are given by 



m l = 



= -3-yHp- -(4-37)^ — — — , 
= -4Hp r , 



i^ 2 +V((f>)+p + p r . 



(26) 

(27) 
(28) 
(29) 

(30) 



where a is the scale factor, H — a/ a is the Hubble parameter, p denotes the energy density 
of barotropic matter, p r denotes de energy density of radiation, <j) denotes the scalar field 
and V((f>) and x{4>) are j respectively, the potential and coupling functions. To maintain the 
analysis as general as possible, we will not specify the functional forms of the potential and 
the coupling function from the beginning. Instead we consider the general hypothesis V{(j)) € 
C 3 , V((f>) > 0, x(</0 G C 3 and x{4>) > 0- We impose they in order to obtain dynamical systems 
of class C 2 . However, to derive some of our results we will relax some of this hypothesis, 
or consider further assumptions (they will be clearly stated when applicable). We consider 
also p > and 0<7<2,7^|. This hypothesis for background matter are the usual. We 
assume 7 7^ | to exclude the possibility that the background matter behaves as radiation. 
The energy momentum tensor for radiation (7 = 4/3) is traceless, so it is automatically 
decoupled from a scalar field non-minimally coupled to dark matter in the Einstein frame. 
Radiation source is added by hand in order to model also the cosmological epoch when 
barotropic matter and radiation coexisted since we want to investigate the possible scaling 
solutions in the radiation regime. We neglect ordinary (uncoupled) barotropic matter. 



3.3 Late time behavior 

The purpose of this section is to formulated a proposition that extent in some way (we are 
considering only flat FRW models) the proposition 1 of [119] , which gives a characterization 
of the future attractor of the system (|2%]) . (|2T|) . pg]) . (|2"9"j) . 

In the following, we study the late time behavior of solutions of (|2"B1) . (|2"T|) . (f2"8")l , (|2"9"j) . 
which are expanding at some initial time, i. e., H(0) > 0. The state vector of the system is 
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(j), (f>, p, p r , H^j . Defining y := </>, we rewrite the autonomous system as 



<t> = y, (31) 

V = - 3Hy -— + 2 {4 ~^ )p ^^> (32) 

P = -37^4(4-3^^), (33) 

Pr = -4#p r , (34) 

# = -UiP+^pr + y 2 ), (35) 



subject to the constraint 



2,H 2 = l -y 2 + V{<i>)+p + p r . (36) 



Remark 1 Using standard arguments of ordinary differential equations theory, follows from 
equations (|33|) and (|34[) i/iai i/ie signs o/ /? and p r , respectively, are invariant. This means 
that if p > and p r > /or some initial time to, then p(t) > 0, and p r (i) > throughout the 
solution. From (|35|) and (|36j) and only if additional conditions are assumed, for example 
V{4>) ^ an d V(4>*) = f or some 4>*> follows that the sign of H is invariant. From ()34p 



and (|35p follows that p r and H decreases. Also, defining e = ^y 2 + V(tf>), follows from 
(|32 ]> -([33 ]) that 

e + p=-3H(y 2 + 1 p). (37) 
Thus, the total energy density contained in the dark sector is decreasing. 

Let us formalize notion of degenerate local minimum introduced in [119 : 

Definition 5 The function V((f>) is said to have a degenerate local minimum at (j)* if 

^'(0),^"(0),...^ (2 ™- 1) 

vanish at 0*, and V^ 2n \(f)^) > 0, for some integer n. 

Proposition 3 Suppose that V((f>) € C 2 (R) satisfies the following conditions^ 
(i) The (possibly empty) set {<f> : V(4>) < 0} is bounded; 
(ii) The (possibly empty) set of critical points ofV(<p) is finite. 
Let 4>* a strict local minimum for V(<j)), possibly degenerate, with non-negative critical value. 
Then p* := [ 0*,y* =0,/?* = 0,p r = 0, H = \ Ylpd \ i s an asymptotically stable equilib- 



3 

rium point for the flow of (|3ip - (|35[) . 
Proof. 

We adapt the demonstration in [119] (for flat FRW cosmologies) to the case where 
radiation is considered. 

3 See assumptions 1 in I119| . 

4 From physical considerations we can neglect radiation for the analysis of the future attractor, but we 
prefer to offer the complete proof. 
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First let us consider the case V ((/>*) > 0. Let V > V(<fr*) be a regular value for V such 
that the connected component of V^ 1 ((— oo, V"]^ that contains 0* is a compact set in M. 
Let us denote this set by A and define as 

9 = {(4>,V, P ,Pr,H) : 4> e A,e + p < V, p > 0, Pr £ [0,W^} , 

where W is a positive constant. We can show that VP is a compact set as follows. 

(i) ^ is a closed set in R 5 ; 

(ii) Vfa*) < V(^) < V", V0e A; 

(iii) iy 2 + V(0*) < \y 2 + V(<j>) + p < V, and therefore y is bounded; 

(iv) p < V — \y 2 — V{(j>) < V — ^ ((/>*) and then p is bounded; 

(v) Finally, from flU, <H 2 <V + W. 

Let C "J" be the connected component of \& containing p*. Following similar argu- 
ments as in |119j it can be proved that ty + is positively invariant with respect to (|3T|) - (f35|) . 
i.e, all the solutions with initial data at ty + remains at ty + for all t > 0. Indeed, let x(£) be 
such a solution and 

t = sup {t > : if (i) >0}elU {+oo}. 

When t < t, equation (157)) imply that e + p decreases monotonically (cf remark[T|). Moreover, 
it can be proved by contradiction that 

4>{i) e A\Jt<t, (38) 

otherwise there would exists some t < t such that V {<)>(£)) > V, but then 

V < V(cf>(t)) < ±y(t) 2 + V{cj>{t)) + P(t) < V, 

a contradiction. Thus, (|38|) holds. But since p r > along the flow (cf remark [1]), it follows 
that 



H{tf > 3 (*) + + >>(*) J ^ ^ir^ ^ ~ip > °- 

We have proved that as long as if remains positive, it is strictly bounded away from 
zero; thus i = +oo, and from this can be deduced that x(i) remains in for all t > 0. 

From all the above satisfies the hypothesis of LaSalle's invariance theorem |126j . If 
we consider the monotonic decreasing functions e + p and p r defined in 'J-), then follows that, 
every solution with initial data at must be such that if (y 2 + jp) —> and Hp r — > as 
t —> +oo. Since if is strictly bounded away from zero in ty + follows that (y, p, p r ) — > (0, 0, 0) 
and H 2 - ^ as t -> +oo. 

Since if is monotone decreasing (cf remark [T]) and it is bounded away from zero it 
must have a limit. This means that V(<f>) also admits a limit. This limit has to be V(<p*); 
otherwise V'((f>) would tend to a positive value and so would the righthand side of (|32[) . a 
contradiction. Therefore the solution approaches the equilibrium point . 

If V(0*) = 0, the above argument can be easily adapted. In this case the set il is 
connected and we choose f2+ to be its subset characterized by the property if > 0. The 
only point in f2 + with if = is exactly the equilibrium point p*, and so if H(t) — > the 
solution is forced to approach the equilibrium since if is monotone; if by contradiction H(t) 
had a strictly positive limit, we could argue as before to find y — > 0, p — > and p r — > and 
so if must necessarily converge to zero. I 
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3.4 Early time behavior 

In order to analyze the initial singularity (and also, the late time behavior) it is convenient 
to normalize the variables, since in the vicinity of an hypothetical initial singularity, the 
physical variables would typically diverge, whereas at late times they commonly vanish 
|161j . In this section we rewrite equations (f^BI . (f2T|). (P251) . as an autonomous system 
defined on a state space by introducing Hubble-normalized variables. These variables satisfy 
an inequality arising from the Friedmann equation ([3"U|) . We analyze the cosmological model 
by investigating the flow of the autonomous system in a phase space by using dynamical 
systems tools. 

3.4.1 Normalized variables and dynamical system 

Let us introduce the following dimensionless variables 



, <ft VP vv ,„ Q , 

a% = <p, <t 2 = —= — , 03 = — p — , 04 = — 7= — , 05 = — 7= — (39) 

Vqh Vzh Vsh Vsh 

and the time coordinate 

dr = 3Hdt. (40) 

We considered the scalar field itself as a dynamical variable. 

Using these new coordinates equations (j26|) - (l29l) recast as an autonomous system satis- 
fying an inequality arising from the Friedmann equation (|30|) . This system is given by 



\°2 (41) 



J j. 1 fw2 _i_ (C\ „ a l dlnX/ ( g i) , (4 ~3 7 ) qf dlnx(txi) 

CT 2 + 7T + — D) (T2 7= ; I 7= • 



4 = ^3 (6.| + 3 7 (4 - 1) + 4.1) - (4-3^^3dlnx(. l) ; (43) 

/ 1 t„ 2 o 2 „ 2\ v 7 ^ dlnl/(cri) 
0^ = -<7 4 (6a 2 2 + 3 7 a| + 4a|) + — a 2 a A j^T"^, (44) 

^ = ^5(6^+37^+4^-4). (45) 

The system (|4T |) -(|45 |) defines a flow in the phase space 

S := I it e I 5 : ^a] = l,aj > 0,j = 3,4,5 ^ . (46) 
i=2 

3.4.2 The topological properties of the phase space 

Let us define the sets So := {0 G £ : 05 = 0} and S + = {<r G £ : 05 > 0}. By construction 
these sets are a partition of E. 

Notations: ||.|| denotes the Euclidean vector norm; D ra := {x G K n : ||x|| < 1} denotes 
the n-dimensional unitary disc; and H™ := {x G M™ : x„ > 0} denotes the n-dimensional 
Euclidean semi-space. 
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Definition 6 ((Topological) manifold with boundary) Let M be a Hausdorff space 
provided with a numerable basis. Let p G M. If there exists a positive number m (possibly 
depending on p), a neighborhood V(p) of p and an homeomorphism h : V(p) — > H m such 
that h(V(p)) C H' m is an open set ofW 71 , then, M is a (topological) manifold with boundary. 

Definition 7 (Boundary) Let M be a manifold with boundary. The boundary of M , dM, 
is defined by dM := {p G M : h{p) £ H m_1 x {0}} . 

Definition 8 (Interior) Let M be a manifold with boundary. The interior of M , IntM, 
is defined by IntM := M \ dM. 

Proposition 4 So is a manifold with boundary. 

Proof. Since So C R 4 is a closed set with respect the usual topology of R 4 , it is 
a Hausdorff space equipped with a numerable basis. The rest of the proof requires the 
construction of a set of local charts. 

Let us define the sets Vj := {a e R 4 : <7j > 0} n So, j = 3, 4. These sets are open with 
respect to the induced topology in So- 

Let be defined the projection maps 

h 3 : V 3 -> H 3 , a -> h 3 (a) = (a 1 ,a 2 ,or 4 ), 

and 

hi : 14 -> H 3 , a -> h 4 (a) = er 2 , <r 3 ). 

These maps satisfy hj(Vj) = K x (H 2 n B 2 ), j = 3,4, which are open sets of H 3 . Their 
inverses are given by 

h^ 1 :IX (H 2 HD 2 ) -> V 3 >1,<72,<7 4 ) -> ^i.tra.^/l-of -o2,o-4^ , 

and 

/14 1 :Rx (H 2 nD 2 )^y 4 , ((Ti,a 2 ,CT 3 ) -> (Vi, <7 2 , <t 3 , -y/l - <r| - a$j . 

It is clear that they are homoemorphisms. 

Observe that {(V3, /13), (V4, /14)} does not cover the sets with 03 = 0-4 = 0. The construc- 
tion is completed by defining the sets 

Vt ■= W G E : cr 2 = ±1, 0-3 = cr 4 = 0} 

which are disjoint copies of R, and thus they are 1-dimensional manifolds. ■ 

By definition (V3, /13), (V4, hi) are topological manifolds. It can be proved that, in fact, 
each one is a topological manifold with boundary. The boundaries are 

dV 3 = {a G V 3 : <74 = 0} , 

and 

dV 4 = {(7 G V4 : (73 = 0} . 

Both are homeomorphic to 1 x (—1, 1). 

Let be defined V x := Vf U Vf. Observe that 

E = (S \ Vx) UVi = Int (S \ Vi) U U SVi) U Vi 

= Int(Eo\^i)U(flEo) 1 U(flEo) a , (47) 

where we have defined (0£ o )i = dV± U Vi = {cr G S : cr 3 = 0} and (<9S ) 2 = <9V 3 U V x = 
{a G S : cr 4 = 0} . 

From the above arguments and expression (|47[) we have: 
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Remark 2 • The interior o/Eo is given by IntT,Q = Int(Y,Q \ V±) which is a 3- dimensional 
manifold (without boundary). 

• The boundary o/Eo is the union of two 2-dimensional topological manifolds with bound- 
ary given by (9£o)i an d (<5Eo) 2 • 

• (t?£o)i an d (9£o)2 share the same boundary V\ which is the union of two disjoint 
copies of M. 

Proposition 5 E + is a topological manifold with boundary. 

Proof. It is easy to check that E + C K 5 is a Hausdorff space equipped with a numerable 
basis. The rest of the proof requires the construction of a set of local charts. 

Let us define the sets Wj := {cr e K 5 : o-j > 0} n E + , j — 3, 4. These sets are open with 
respect to the induced topology in £ + . 

Let us define the maps 

g 3 : W3 -> H 4 , cr -> 5(3(0-) = (01,02,05,04) , 

and 

34 : W 4 -> H 4 , o -> 54(0-) = (0-1,02,05,03) . 

These maps satisfy g 3 (Wj) = K x (H 3 nB 3 ), j = 3, 4 which are open set of H 4 . Their inverses 
are given by 



g 3 1 :lx (H 3 n D 3 ) ->■ W3, (01, o 2 , o 5 , o 4 ) foi, o 2 , y'l - of - a\ - of, 04, 05 

and 



g 4 1 : K x (H 3 nD 3 ) -> W4, (01,02,05,03) -> (^01,02,03, yi - o| - of - of, 05 

It is clear that they are homeomorphism. 

Observe that {(^3,33), (W4, 5/4)} do not cover the sets with 03 = 04 = 0. The construc- 
tion is completed by defining the set W\ := {o 6 £ + : o^ + a\ = l} . Using the projection 

map (01,02,05) ^\ (01,02) it is easily proved that Wi is homeomorphic to the open set 
E x (-1,1). ■ 

We have that W3 is a manifold with boundary. Its boundary is the set dW% := {o e E : 
03 > 0, 04 = 0, 05 > 0} and its interior is the set {o G S : 03 > 0, 04 > 0, 05 > 0}. Also, W4 
is a manifold with boundary. Its boundary is the set dWi := {o € E : 03 = 0, 04 > 0, 05 > 0} 
with same interior as W3 and W\ is a manifold without boundary which is homeomorphic 
to R X (-1, 1). 

Let us define the sets 

(9S+) X := <91U 4 UVUi 

and 

(<9E+) 2 := dW 3 U 1Ui. 
Following the same arguments yielding to (|47|) we get the formula 

£+ = Int (£+ \ Wi) U (9S+)! U (5E+) 2 . (48) 

Using the above arguments and relation ({48} we have the following 

Remark 3 • The interior of E + is given by IntT, + = Int (E + \ W\) which is a 4~ 
dimensional manifold (without boundary). 
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• The boundary of £+ is the union of two 3- dimensional topological manifolds with 
boundary given by (c?£+) 1 and (d£+) 2 . 

• (9E + ) 1 and (<9E + ) 2 share the same boundary W\ which is a 2-dimensional manifold 
without boundary homeomorphic to R x (— 1, 1). 

3.4.3 Monotonic functions 

The construction of monotonic functions in the state space is an important tool in any phase 
space analysis. The existence of such functions can rule out periodic orbits, homoclinic 
orbits, and other complex behavior in invariant sets. If so, the dynamics is dominated 
by critical points (and possibly, heteroclinic orbits joining it). Additionally, some global 
results can be obtained. In table [1] are displayed several monotonic functions for the flow of 
flH])-(|45]) in E. 



Table 1: Monotonic functions for the flow of ([4"T ]) -([4"5 |) in (|4"5|) . 



z 


dZ/dr 


Invariant set 


Restrictions 3, 








a 3 > 0, 04 > 


7 6(0,§)U(§,2) 






-(!-7)^ 2 


(T 3 > 0, f7 5 > 


7 e(o,|)u(|,2) 








04 > 0, 05 > 


none 




-§Z 4 


C2 7^ 0,0-3 = 0-4 = 0, 05 > 


none 



a We assume the general conditions %, V € C 3 , x(ai) > 0, V^(cti) > 0. 



Remark 4 From the definition of Z\ it follows that it is a monotonic decreasing func- 
tion for the flow of (|41[) - (|45|) restricted to the invariant set 03 > 0,04 > 0. Applying the 
Monotonicity Principle (theorem^ it follows that the past attractor the flow of (|41[) - (|45[) 
restricted to the invariant set 03 > 0, 04 > is contained in the set where 04 = and the 
future attractor in contained in the set where 173 = 0. 

Remark 5 Using the same argument follows from the definition of Zi that the past attractor 
of the flow of (|41[) - (|45[) restricted to the invariant set 03 > 0,05 > is contained in the 
invariant set where 03 = and the future asymptotic attractor is contained in the invariant 
set 05 = provided 7 < If 7 > | the asymptotic behavior is the reverse of the previously 
described. 

Remark 6 From the definition of Z3 follows that the past attractor of the flow of (|41[) - (|45[) 

restricted to the invariant set 04 > 0, 05 > is contained in the invariant set where 04 = 
and the future asymptotic attractor is contained in the invariant set 05 = 0. 

Remark 7 From the definition of ' Z4 follows that the past attractor of the flow of (|41[) - (|45[) 

restricted to the invariant set 02 ^ 0, 03 = 04 = 0, 05 > is contained in the invariant set 
where 02 = (i.e., where 05 = 1) and the future asymptotic attractor is contained in the 
invariant set 02 = ±1. 
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3.4.4 Critical points with <f> bounded 



Let us make a preliminary analysis of the linear stability of the critical points of the flow 
of (|4"TT) - (|4"5)) defined in S. It is a classic result that the linear stability of the critical points 
does not change under homeomorphisms. 

Since E is a 4-dimensional manifold (with boundary) we will consider the projection of 
E in a real 4-dimensional manifold (with boundary). 

Let be defined the projection map 



7T : E -> n 

(<Tl, C72 , 03, CT4, CT5) — > (<7l, <72, 03, CT5) 



(49) 



where 



O := {(i £ R 4 : cr 2 2 + + (T5 < l,aj >0,j = 3,5} . (50) 
The flow of (|4"Tj) - (j4"5|) defined on E is topologically equivalent (under n) to the flow of 



n ~< = \l -n> (51) 

, ^3, l/o 2, 4 2 r n (l-^-^3 2 "^)dln^i) , (4-3 7 )a 3 2 dlnx(^) 
- 2 = - 2 + g(37a 3 +4a 5 -6)a 2 __ + __, 

(52) 
(53) 

(54) 



1 



4 = g ^3 (6a 2 2 + 3 7 (a 2 - l) + 4a 2 ) - 



(4- 37) CT2CT3 dlnx(cri) 



diTi 



a' = ^a 5 (6a 2 + 3 7 a 3 2 + 4a 2 - 4) , 



1 

defined in 



Table 2: Location of the critical points of the flow of (f5"T j) - (f54"l) defined in f2. 



Label 




02 


CT3 


05 


Qi 


(Tic ■ X'(°"Lc) = 





1 





Q2 


vie : V'(<t 1c ) = 











Qs 


<Ti c e K 








1 



We are assuming V(ai c ) ^ 0, x( tJ ic) 7^ and 7 € (0, 1) U (1,2) 



The system ([5T|) - (f54|) defined in fi admits three classes of critical points located at O, 
denoted by Qj, j — 1,2,3. In table [2] are displayed the coordinates of such points. The 
dynamics near the fixed points (and its stability properties) is dictated by the signs of 
the real parts of the eigenvalues of the Jacobian matrix evaluated at each critical point as 
follows: 

1. The eigenvalues of the linearization around Qi are — A2, 7, Ai ± ^A 2 + A2 \ 

where Ai = (-2 + 7 )/4 < 0, and A 2 = (4 - 37) /6. Then, the local stability of critical 
point P\ is as follows: 
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(a) If — Al ^^ lc ' ) < x"( a ic) < and < 7 < | there exist a 3-dimensional stable 
manifold and a 1-dimensional unstable manifold of Q\. 

(b) If < 7 < § and x"(<?ic) > or § < 7 < 2 and < x"(cr lc ) < - A 'ffi le) , there 
exist a 2-dimensional stable manifold and a 2-dimensional unstable manifold of 

Qi. 

(c) If < 7 < I and x"i cr ic) = there exist a 2-dimensional stable manifold, a 

1- dimensional unstable manifold and a 1-dimcnsional center manifold of Q\. 

(d) If I < 7 < 2 and x"{ a ic) = there exist a 1-dimensional stable manifold, a 

2- dimensional unstable manifold and a 1-dimensional center manifold of Q\. 

2. The eigenvalues of the linearization around Q 2 are: — ^ , — ^ , — | ± — | \/(^ 1C )' > ■ 
Then, the local stability of critical point P 2 is as follows: 

(a) If V"(cti c ) < there exists a 3-dimensional stable manifold and a 1-dimensional 
unstable manifold of Q 2 - 

(b) If V" ((Tie) = the eigenvalues are — 1, — |, 0, — ^. Then there exists a 1-dimensional 
center manifold tangent to the 01-axis and a 3-dimensional stable manifold tan- 
gent to the 3-dimensional surface 

< (<7i,<7 2 , 0-3,0-5) € M 4 : 01 = 



(c) If < V"(<7i c ) < I V(oic) the stable manifold of Q2 is 4-dimensional (the critical 
point is asymptotically stable and it is an stable node). If V"{<Ji c ) > fV(oi c ) 
the stable manifold of Q 2 is 4-dimensional (the critical point is asymptotically 
stable and it is an stable focus). 

3. The eigenvalues of the linearization around Q3 are 3,-3, 0, g(4 — 37). Then, the local 
stability of critical point Q3 is as follows: 

(a) There exists a 1-dimensional center manifold tangent to the 01-axis 

(b) If 7 < |, Qz has a 2-dimensional unstable manifold and a 1-dimensional stable 
manifold tangent to the line 

l (01,02,03, 05) e K 4 : 01 = -^o 2 , o 3 = o 5 = j . 

(c) If 7 > |, Q3 has a 1-dimcnsional unstable manifold and a 2-dimensional stable 
manifold tangent to the 2-dimensional subspace 




(oi,o 2 ,o 3 ,o 5 ) e K 4 : 01 = |cr 2 , o 5 = 

Let us comment on the physical interpretation of the critical points listed above. The 
critical points Qi represent matter-dominated cosmological solutions with infinite curvature 
(H — > +00). Since x'(oic) = (i.e., o\ c is an stationary point of the coupling function) 
and x'(oi c ) 7^ they are solutions with minimally coupled scalar field and negligible kinetic 
energy. The potential function does not influence its dynamical character. 

The critical points Q 2 represent a de Sitter cosmological solution. When one of these 
critical points is approached, the energy density of dark matter and the kinetic energy of 



22 



the scalar field go to zero. In this case the potential energy of the scalar field becomes 
dominant. Hence, the universe would be expanding forever in a de Sitter phase. 

The critical points Q3 represent radiation-dominated cosmological solutions. These so- 
lutions are very important during the radiation era. 

3.4.5 Normal form expansion for (|5Tj) - (f54|) around Q2. 

In this section we present normal form expansions for (|51|) -(|54 j) around Qi- We exclude 
from the normal forms analysis the case where V'(ai c ) = and V"(ai c ) > (i.e., V((f>) 
has a local minimum at a\ c ) since from the linear analysis (see point 2c in subsection 13.4.4( 1 
follows the asymptotic stability of Q2- 

Let V(4>) and x(0) De smooth (C°°) functions. The barotropic index of the fluid satisfies 
0<7<2 )7 ^|. 

Case 1. Let a lc such that V'{a lc ) = V"(a lc ) = and V^{a lc ) ^ 0.@ 

Let us consider the coordinate transformation {a\, 02, 03, 05) — > (cri + a\ Cl 02, 03, 0-5). 

The Taylor expansion up to third order of the of the system arising from (|5I[) - (f54"|) under 

such coordinate transformation around the origin reads 



3^2 + 0(3), 



(T 2 = -CT 2 - 



o\ V^(a lc ) , (4-3 7 )a| X / (a lc ) 



2^6 V(a lc ) 2V6 x(o-ic) 

7 (4 - 37)^2^3 x'(°"ic) 

|<7 B + 0(3). 



2V6 X(-lc) 



0(3), 



0(3), 



(55) 



1. Let be 7 ^ 1. 

Applying the successive coordinate transformations 



/ -1 ^ 




02 




0-3 




V J 


V 



0-3 



-1 



01 

05 

(T 2 



/ 



(56) 



-3 
V <T 5 J 



( n (3 7 -4) X '( CTle yg (2^-3<t 3 2 )v< 3 )( CT1c ) \ 
1 " t " 2v / 6(7-l)x(<Tie) 6V6^(<Tie) 



03 



V 



(72 

(3 7 -4) x , (<T le )g-g cri(o- 1 -276CT 3 )v' 3 >(<Ti c ) 
67x(^ic) + 18V(ffic) 

_ I (4-37)q-l (T5X'(o , lc) 

175 + 2v^6 x ( CTlc ) 



The system ([55]) becomes 



(57) 



5 In this case a\ c is not an extremum point of V(<j>) but an inflection point. 



23 



01 



03 



y0)( 



(Tic 



3V(a lc 
-^ + 0(3), 



0(3), 



6^(ai c ) 

2 



■0(3), 



0(3). 



(58) 



Neglecting the terms of third order in (|58|) the exact solution of the approximated 
system passing through (crio, c 2 o, cr 30 , (T50) at t = r is given by 

' (T-T )cr 30 VW(a lc y 2 



(71 (t) 


= e- (r - ro Vio ( 


(72 (t) 




cr 3 (r) 


= <7 3 ^1 + 


ct 5 (t) 





6F(a lc 



From solution (15^1) it is easy to check that 



lim (<j 1 (T),a 2 (T),a 3 (T),a 4 (T)) = (0,0,0,0) 



(59) 
(60) 



at the rate (1 + 0(r)) 1 . This means that, under the conditions discussed here, the 
origin of coordinates is asymptotically stable for the flow of the original system ([55]) . 

2. Let be 7 = 1. Applying the successive coordinate transformations 



( ax \ 

0"2 
0-3 

V / 



05 



2 „ 
3CTI 



V 



01 
03 
0~2 



/ 



CT3 
V CT 5 / 



Cl 



V 

The system (1531) becomes 



°5 + To 



(2^-3 CT5 2 )y' 3 )(er le ) 
6v / 6l / (<Ti c ) 

I <7lT3X'(glc) 

CT3+ 2VSx(oia) 

<T 1 (<T 1 -2V6<T S )V (3) (Via) _ 3q|xVle) 
V(cri c ) X(cic) 



0"l = O"! 



0-5 



F (3)( 



CTic 



3F(ai c 



- 1 



XVle)gf 

2\/6x(aic) 



0(3), 



0"2 



2ff2 

3 



0(3), 
?( 

gfV (3 Vic 

6V(<7lc) 



-? + O0), 



0(3). 



(61) 



(62) 



(63) 
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Neglecting the terms of third order in (1551) the exact solution of the approximated 
system passing through (<7io, 020 j 030, 050) at r = To is given by 

2 



a 1 (r) = e-( r -^<T 1Cl ( I 



(T-T )a 50 V^(a lc ) 



(t - t )x' (picMo 



<Ji (r) 

5 (r) = ( 1 



01c 



e-t^-^o-ao, 



6F(a lc ) 



From solution (|64|) it is easy to check that 



lim (cri(r),CT2(T),CT3(T),CT 5 (r)) = (0,0,0,0) 



(64) 



(65) 



at the rate (1 + 0(t)) 1 . This means that, under the conditions discussed here, the 
origin of coordinates is asymptotically stable for the flow of the original system ([55]) . 

Case 2. Let the function V(4>) to have a degenerate local minimum at u\ c (for n = 2). 

Let us consider the coordinate transformation (01, 02, 03, 175) — > (ci + a\ c , 02, 03, 05)- 
The Taylor expansion up to fourth order of the of the system arising from (f5"Tj) - (f54"|) under 
such coordinate transformation around the origin reads 



0i = \/ 302 + 0(4), 



VW(vic)al | (3 7 = j)a|xVic)V , (4 - 3 7 )0 3 2 X>ic)0i 
6V6F( ( n c ) 



02 - 02 + g02 (370-3 + 40^ 



2 v / 6x(0ic) 2 2v / 6x(0i c ) 
~ (4-3 7 )a 3 2 X '(a lc ) 



(4 - 37)ctict 2 03X'(01c) 5 



2V6x( 



01c 



2V6 X (0ic) ' ° (4) ' 
(3 7 -4)cr 2 er 3 x'(0i c ) 703 
2 



6 



a 3 (6o% + 3 7 a| + 4cr|) 



2\/6x(0i c ) 
(3 7 - 4)gia 2 03X"(0ic) 
2V6x(0-ic) 



•0(4), 



1 



05 = g05 (60 2 + 37CT 3 + 4(T 5 

1. Let be 7 7^ 1. Applying the successive coordinate transformations 



^ + 0(4). 



(66) 



fa, \ 

02 
03 

V °* J 



02 
03 
V 05 / 



^ 03 ~ 

V 



2 „ 
3CTI 



01 



03 



01 

05 
02 



V 05 



(37-4)x'(<ri<>g 
2V6( 7 -l) X ( ' lc ) 
02 

_ (3i-±)x' (iic)°l 

67X(°'lc) 
(4-37)crio-5x'(cic) 
2V6x(^ic) 



(67) 



(68) 
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V ^ J 
( 



/ 2 , 2 , 2\ (2aj i -6v / 6cr 3 ^+3v / 6crij)V (4) (<n e ) \ 
1°"! + °2 + CT 5j 108VTW) 



<J 1 [-2a l f 3 —^ 

-I 1 

S 3 - (a? + a\ + o\) 



<Ti(2V^CTi-27cr 3 (Ti+27v / 6(T3)V <4) (<Ti c 
486V(cri c ) 

(37-4)cri (cti-v / 6ct 3 )(T5x"( '1c) 



(69) 



where 

„ _ (37-4)(3((7- 3 )7+2)cr 1 +V6 7 a 3 ) x , ( C r le ) 2 t Tg _ (3 7 -4) (2(7- lVi-y^^x" (eri e ) 2 
" 1 ~~ 12(7-1) 7X (<t 1c )^ 12( 7 -1)7x(<t 1c ) 

„ _ (37-4)(- ^f" -^)x> lc ) 2 ^ 2 (37-4)(^-^)x"(^c) g5 2 
"2 - 12 x (<xi c ) 2 18 X (<rlc) 

and 

_ (37-4)<T 5 (8g- 2 + 7 ((7-l)<Ti((37-8)<Tl+476 C r3)-6g- 2 ))x , (cri e ) 
~ 48( 7 -l)7x(o-i c ) 2 



-3 



the system l|6"6")) becomes 



■^ + 0(4), 



V(a lc 
O 

18V(a lc ) 
-2^ + 0(4). 



■0(4), 



■0(4), 



2. Let be 7 = 1. Applying the successive coordinate transformations 



( CT1 ^ 




(0* 






— ► 




CTl 


^3 






^3 


V ^ ) 




\ 


0-2 / 





^1 










02 










^3 








V 


CT5 


J 







0-3 



0"! 
02 

q"ig"3x'(^ic) 

2V6x(tic) 
6x(o"ic) 



(70) 



(71) 



(72) 
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(*1 \ 

CT"3 

V °* J 



/ 



7? 1^- /,t2 



(2 g -;-6v / 6g 5 ^+3v / 6cri)V (4) (<n e ) 
108V(cti c ) 



7^2 (erf 



77 i^- /^2 

=5 - 2 CT 3 ^1 

o-i(28crf +24cr|+21o-: 



cr. 



103 (fi - V6o-5 )x" ( CT lc ) 
12x(o-io) 



o-i (2 -27ct 5 cri +27V6cr 2 ) V w (a lc ) 



where H4 = 

(V6cri-4 f 75) x , ( g -i c ) 2 a -| 
24 X (^lc) 



42V6 

1 RvCrr, ' 



4x(^lc) 2 ^ 6x(<Tlc) 

( \/6(Ti -6(T 5 ) x" (trie ) o- 2 
36 X (o-ic) 



486V(cti c ) 

cr 3 (5cr 2 — 4\/6o"5 cri +2(T 2 ) x' ("'lc 
~~ 48x(cti c )2 

the system (frJrJ)) becomes 



(73) 



and S« 



f x>ic) x'(^c) 2 y T 2 , ^ 5 2 ^ (4) (^c) 1S (0(h 

-^+0(4), 
-y +0(4), 



^ (4 Vic) 
W(<r lc ) 



0(4). (74) 



Having at hand the normal forms expansions ((70)) and (l74l) we can extract the following 
information. Let ^(cti) to have a strict degenerate local minimum at (j\ c (for n = 2) with 
V(°ic) > 0. Let 7^1. Neglecting the terms of fourth order in (|TP1) the exact solution of 
the approximated system passing through (ctio, 020, "so, CT 5o) at t = To is given by 



CT 3(^) = 77777 \ 1" 1 O"30, 



ai(r) = l m^) + l ) ' n] '" 

(72 (t) = e-t( T - T °)(720, 

( r _- ro ) F (4)( gl e)(7 | 

9K(ai c ) 

a 5 (r) -e-* (r - ro V 50 . (75) 
From solution (fT5"j) it is easy to check that 

lim (a 1 (T),a 2 (T),a a (T),a 5 ( T )) = (0,0,0,0) (76) 

at the rate (1 + O(t))- 1 / 2 . 

Let 7 = 1. Neglecting the terms of third order in (f74| the exact solution of the approxi- 
mated system passing through (ctio, 020, 030, (750) at r = tq is given by 



27 



( 



9V(a lc ) 




V 2 v / 6x(ai c ) 2 
a 2 (r) = e-i( T - r °V 20 , 
o-a (r) = e-^- TO V 3 o, 



/ (t - T)a 5 ox , (gic) 2 g| ( 



+ 010 + 



(r-r ) ( 7 1 o ( T 5 2 o y( 4 )(a le ) 
9V{a lc ) 



) 




50 



+ 1 



) 



2 



^50- 



(77) 



From solution ([77)) it is easy to check that 



Urn (<ri (r), <r 2 (r), (t),ct 5 (t)) = (0,0,0,0) 



(78) 



(l + oW)-'/ 2 . 



Hence, for all 7 G (0,4/3) U (4/3,2), the origin of coordinates is asymptotically stable 
for the flow of the original system ([SrJ| . Since the remainder of the Taylor expansion (|rIS|) is 
0(|er| 4 ) as cr — > 0, uniformly in To < t < 00, follows the asymptotic stability of the origin. 
This analysis is in agreement with the result of proposition [3] which states that if V(a±) 
have a strict degenerate local minimum at o\ c with V{a\ c ) > 0, then Q2 := (cti c , 0,0,0) is 
asymptotically stable. 

3.4.6 A generalization of theorem 3.2 of [lj p. 8. 

Theorem 3.2 of [T| p. 8, states, essentially, that if the potential and the coupling function 
are sufficiently smooth functions, then for almost all the points p lying in T,q, then the scalar 
field diverges through the past orbit of p. This result can generalized as follows. 

Generalization 

Theorem 8 Assume that x(0) an d V((f>) are positive functions of class C 3 , such that x has 
at most a finite number of stationary points and does not tend to zero in any compact set of 
R. Let 7 £ (0, |) U (|, 2) and let p be a point in T, + . Let 0~(jp) be the past orbit of p under 
the flow of (|41[) - (|4"5l) restricted to S + . Then, <f> is always unbounded on 0~{p) for almost 
all p. 

Proof of theorem \S\ In order to prove the theorem it is sufficient to consider interior 
points of E+. Also, in order to apply results concerning future attractors and w-limit sets 
we perform the time reversal r — > — r. Thus we get the system 
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2 <7 2 (79) 

6 v ' 3 5 ' V6 do-i 2^ doi 

' 1 /« 2 , o / 2 ,\ , , 2\, (4-37) g 2 ^3 dlnx((Ti) 
03 = -gff3 (6a 2 + 3 7 (o- 3 - 1) + 4a 5 ) + ^= — , (81) 

^4 = - g^4 (6<r 2 2 + 3 7< 7| + 4of ) - ^-0204 v u , (82) 
05 = - ^5 (6<t| + 3 7 a| + 4tr| - 4) . (83) 

where the prime denotes now derivative with respect to — r. 

Let po := (<Tio, (T20, 030, 040, "so) G IntS+ such that there exist a real positive number K 
with |<7i| < K for all p := (01, <j 2 , 03, (T4, og) <E + (po), where + (po) denotes the positive 
(future) orbit for the flow of ([791)-®. Then, for all p G O+(p ) we have 

-1 < a 2 < 1, < cr 3 < 1, < cr 4 < 1, < cr g < 1. 

Hence + (po) is contained in a compact set of (the closure of) £+. 

Since + (po) is a positive invariant set, then using proposition [I] we ensure the existence 
of a non empty, closed, connected and invariant w-limit of po denoted by w(po)- 

First we demonstrate by contradiction that that 03 and 174 cannot be simultaneously 
zero at w(po)- Suppose that w(po) is contained in the set where 03 = 04 = 0. Let us define 
the function Mi = Z^ 1 (see table [T] for the definition of Z 4 ) defined in the invariant set 
&3 = 04 = 0, < 05 < 1. From the definition of Mi and applying the Monotonicity principle 
(theorem [T]) follows that the future asymptotic attractor of the flow of (|79 p -(|83 p restricted 
to the invariant set (73 = er 4 = is contained in the invariant set a% = ±1. Thus, from (|79[) 
follows that <f) — > =Foo as w(po) is approached, a contradiction. 

Second, let be defined in 7ntE+ the function M 2 = Z^ 1 (see tabled] for the definition of 
Zi). The derivative of M 2 along any orbit of (f79" ]) -([83 ]) is given by M 2 = -7M2. Then M 2 is 
a C 3 monotonic decreasing function for the flow taking values in the interval (0, +00). Since 
03 and 04 cannot tend to zero simultaneously in uj(p) for p G + (po), then the function 
M 2 tends asymptotically to a well defined limit. By construction M 2 (j>) — ^ if and only if 
p — » q with g satisfying cr 4 = (we are using here the condition that x{4>) does not tend to 
zero in any compact of K) and M 2 (p) — > +00 if and only if p — > q with g satisfying (T3 = 0. 
Thus, applying the Monotonicity principle (theorem [1} follows that 

w(po) C {p G S+ : |cri| < if, 0-3 > 0,cr 4 = 0} = Si. 

Let go G u)(po). By the invariance of the w-limit set follows that a; (go) — w (Po)- 

Observe that 33 (Si) = {0" G W3 : \<Ji \ < K, ct 4 = 0} where 33 is defined in the proof of 

proposition 

Let us define the projection map 

g : (01, 02, 03, 05) -t (01, 02, 05) 
and let ctq = .g o 33(go) then the flow of ([79"]) - ([53")) in a neighborhood of g contained in Si, 
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is topologically equivalent to the flow of 



CT 1 = -V {J* 7 *, 



' 1 h v\(n ^ (4-3 7 )dln X (a 1 ) ^ (4 - 3 7 ) / >/6dln X (^) 
^2 = 7t(1- o-a) ( 2 - 7)^2 7= t- t <7 2 



2 v 27 V \/6 do-x / 6 I 2 dax y ' 

4 = -^5 (3(2 - 7 )of - (4 - 3 7 )(1 - a 5 2 )) , (84) 

in a neighborhood of ao contained in 

S = {(crx,cr 2 ,cr 5 ) : -K < <r t < K,al + al < l,a 5 > 0} . 

Since the vector field is C 2 we can extent the flow of (|84|) to the closure of S (denoted by 
S). 

Let us investigate all possible compact, non empty, and connected invariant sets of (|89|) - 
([50)1 located in the closure of S (these ones can be candidates to the w- limit w(<7o)). 
Let us consider two cases: 

i) < 7 < |. Let be defined in S, the function 

CT 5 

The derivative of M3 through an arbitrary orbit of (1841) is given by 

M3 = -2 Q - 7) M a . 

Then Af3 is a C 3 monotonic decreasing function for the flow taking values in the 
interval (0, +00). By construction M^(p) — ¥ if and only if p — ¥ q with q satisfying 
0"2 + erf = 1 (since x does tends to zero in [-K, K]) and M^(jp) — > +00 if and only if 
p — > q with q satisfying 05 = 0. Thus, applying the Monotonicity principle (theorem 
[lj follows that 

w(ff )c{(Te5\ 5:(7 2 2 + (T 5 2 = l}. 
Let go € a; (00). By the invariance of the oj- limit follows that w(cro) = w(qo)- 
Let us define the projection map 

g' : (a-x,cr 2) cr 5 ) (o-x,c 2 ) 

and let cr = g'(qo) then the flow of (|M|) in a neighborhood of go contained in S 1 , is 
topologically equivalent to the flow of 

a> 2 = -(7 2 (1 - <7 2 ) (1 + a 2 ) . (86) 

in a neighborhood of <7 (contained in S" := (—K,K) x (—1, 1)). 
Let be defined in 5" the function 



M 4 (ct) 



l-al 
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which satisfies M{ = — I-M4 along an arbitrary orbit of Thus M4 is a C 3 

monotonic decreasing function in S' . Applying the monotonicity principle (JXJ) follows 
that 

w(^)c{d6 5'\5':(7 2 2 =l}. 

Thus to(cr' ) is contained in one of the invariant sets of (|89p - (|9"0"|) given by a 2 = ±1 but 
this would imply the divergence of <f>. A contradiction. 



ii) I < 7 < 2. Let be defined in S, the function 

M 5 (a) = —^— -j. (ST. 



^ 5 2 x(^i) 37 - 4 



The derivative of M5 through an arbitrary orbit of is given by 

M^ = -2 f 7 _lW. 

Then M5 is a C 3 monotonic decreasing function for the flow taking values in the 
interval (0, +00). By construction M§(p) —¥ if and only if p — ¥ q with q satisfying 
05 = (since \ does not tends to zero in [— K, K]) and M$(p) —> +00 if and only if 
p — > q with q satisfying a\ + c| = 1 . Applying the Monotonicity principle (theorem [T]) 
follows that 

U)(a )c{a€S\S:a 5 =0}. 
Let (jo G w((7o). By the invariance of the w-limit follows that w(cto) = w(go)- 
Let us define the projection map 

h : (ai,a 2 ,<r 5 ) (ai,a 2 )- (88) 

Let a' Q — h(q ). Then, then the flow of in a neighborhood of q$ contained in S, is 
topologically equivalent to the flow of 



3^2, (89) 

rr> l h rr*\((1 „W (4 ~ 3 7 ) d 111 xfo) ^ 

^ = 5(1-^)^(2-7)^--^ j ' (90) 

in a neighborhood of <7 (contained in S"). 

Let us investigate the possible compact invariant sets of (|89I) - ([9H|) located in the closure 
of 5' which can be candidates to the w-limit u;(er ). 

First oj(<j' ) cannot be contained in the invariant sets of (159"l) - (l9T)l) given by a 2 = ±1 
because this would imply the divergence of </>, a contradiction. Second, combining the 
results of the Poincare-Bendixon Theorem (theorem^ and Dulac's criterion (theorem 
13]) with = (1 — cr|) _1 follows that the only possible compact invariant sets are 

the critical points with a\ bounded (or heteroclinic orbits joining such critical points) . 

Let us consider x( a i) other than exponential. In this case the system (|59")) - (f9"0")) 
admits a (possibly empty) family of critical points 

Q := {(qi, 0) e [-K, K] x [-1, 1] : X '(qi) = 0} . 



3 As we will see next in section [4. II the following analysis applies also to exponential coupling functions. 
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If Q = 0, i.e., x'iQi) f° r an l<7i| < A, then the future orbit + (<r n ) tends to a 
point with a\ = ±1. From this follows that tf> is unbounded (a contradiction) and the 
proof is done. 

Let us assume that Q ^ 0. Let q € Q. The eigenvalues of the Jacobian matrix 
g£(?),i,j = l,2are fi ±=A 1 ± ^A? + A 2 ^, where A x = V > °> A ? = ^ 
Hence, at least one of its associated eigenvalues has positive real part. Let be defined 
the sets Q ± = {q € Q : ±x"(<?) > 0} and Q° = {q e Q : x(g) = 0}. At least one of 
these sets is not empty. Let be define A = {p e [—if, A"] x [— 1, 1] : linv^oo <& T (p) = (?}. 
There are the following cases 

• q € Q + , | < 7 < 2, then is 2-dimensional implying A = 0. 

• 5 € Q~, | < 7 < 2, then £"((?) is 1-dimensional and f s (<z) is 1-dimensional. 
Then A C A, leb(N) = 0. 

• g G Q , then £ c (q) is 1-dimensional and £ u (q) is 1-dimensional in such way that 
A c £ c (g), leb{S c (q)) = 0. 

Therefore, all solutions future asymptotic to q (and then with <p bounded towards the 
future) must lie on an stable manifold or center manifold of dimension r < 2, and then 
contained in a subset of [—A, A] x [—1,1] with zero Lebesgue measure. Since there 
are at most a finite number of such q the result of the theorem follows. I 

Theorem [8] allow us to conclude that in order to investigate the generic asymptotic 
behavior of the system (|4T j) - (j45j) restricted to E + it is sufficient to study the region where 
<f> = ±00. 

3.4.7 Analysis in the limit <f> — > oo. 

In this section we will investigate the flow as <f) —¥ oo following the nomenclature and 
formalism introduced in [113] (see also |162j ). Analogous results hold as cj) —> — oo. 

Definition 9 (Function well-behaved at infinity [H3]) Let V : M -> R be a C 2 non- 
negative function. Let there exist some 4>q > for which V{4>) > for all 4> > 0o an d some 
number N such that the function Wy ■ [4>o, oo) — > M, 

satisfies 

lim W v {4>) = 0. (91) 
Then we say that V is Well Behaved at Infinity (WBI) of exponential order N . 

It is important to point out that N may be 0, or even negative. Indeed the class of 
WBI functions of order is of particular interest, containing all non-negative polynomials 
as remarked in |113j . 

In order to classify the smoothness of WBI functions at infinity it is introduced the 
definition 

7 We denote the stable, unstable and center manifolds of q by £" (q), £ u (q) and £ c (q) respectively. By 
leb(A) we denote de Lebesgue measure of i C I 2 . 
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Table 3: Simple examples of WBI behavior at large (p. n and A are arbitrary constants. 



Adapted from |113j . 



V((p) 


Wv((p) 


V — I VP) 


WvvP) 


/ VP) 






<j)~ 2 


ntp 2 









r 1 













Xtp 2 




(A + (<f> - B) 2 ) e-"* 


2(<f>-B) 


0~ 2" 


2 ¥ > 2 (Bip 2 -l) 












2A 2 




2A 2 




l-e A2 « 




A 2 

1-e v 




Oln^)- 1 


(ln^)" 1 


l 


2 


cp 2 ln0 


+ (01n 


(ln^)- 1 


(2 + t^e^ 


2 

—ipe v 



Definition 10 Let be some coordinate transformation tp — f{(p) mapping a neighborhood of 
infinity to a neighborhood of the origin. If g is a function of (p, g is the function of tp whose 
domain is the range of f plus the origin, which takes the values; 



9VP) = 



lim^oo g{cp) , (p = 



Definition 11 (Class k WBI functions [113J) A C k function V is class k WBI if it is 
WBI and if there exists (f>o > and a coordinate transformation tp = f((f>) which maps the 
interval [<po, oo) onto (0, e], where e ~ f((f>o) and lim^co f — 0, with the following additional 
properties: 

i) f is C k+1 and strictly decreasing. 

ii) the functions WvVP) and f'if) are C k on the closed interval [O.e]. 

dtp dtp 

We designate the set of all class k WBI functions £ £ . In table [3] are displayed simple 
examples of WBI behavior at large (p. 

By assuming that V, x € £ 3. , with exponential orders N and M respectively, we can define 
a dynamical system well suited to investigate the dynamics near the initial singularity. We 
will investigate the critical points therein. Particularly those representing scaling solutions 
and those associated with the initial singularity. 

Let S e = {(eri , 02, o~3, 04, 05) G S : o"i > e where e is any positive constant which is 
chosen sufficiently small so as to avoid any points where V or \ — 0, thereby ensuring that 
Wv(<p) and W x (tp) are well-defined, d 

Let be defined the projection map 



tv 1 : S e — » S1 E 

(cri, 02, 0-3, 0-4,0-5) (0-1,0-2,0-4,0-5) (92) 

where 

fi e := {o <= M 4 : 01 > e _1 ,of +<7| +of < l,Oj > 0,j =4,5}. (93) 



8 See ll0l for the definition of functions with bar. 
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Let be defined in f2 e the coordinate transformation (<7i, 02, 04, 175) — > (ip, 02, 05) 
where f(ci) tends to zero as a\ tends to +00 and has been chosen so that the conditions 
i)-iii) of definition [TT1 are satisfied with k = 2. 

The flow of (|4"Tj) - (|4"5|) defined on E £ is topologically equivalent (under / 07^) to the flow 
of the 4-dimensional dynamical system 



V' = y|/V 2 , (94) 
2a 5 2 \ {W v + N)al 

^ / 7! + 

+ [ 27! + -2-J(l-* a -"4-*6) (95) 

(7 4 = ia 4 (V6 (Wy + iV) a 2 + 3(2 - 7 )of + 3 7 (1 - of) + (4 - 3 7 )a|) , (96) 
a' 5 = ±* 6 (3(2 - 7 )af - 3 7 a 4 2 - (4 - 3 7 )(1 - a 2 )) , (97) 

defined in the phase space 

a = {>, (72,(74,(75) £l 4 :0<^< !{^ V )A +al +al < 1,<74 > 0, as > 0}. (98) 

It can be easily proved that defines a manifold with boundary of dimension 4. Its 
boundary, d^, is the union of the sets {p € f2 e : (/? = 0}, {p € f2 e : ip = /(e -1 )}, {p £ Sl ( : 
o"4 = 0}, {p € Q e : (75 = 0} with the unitary 3-sphere. 



Critical points of the flow of ([MJl-flHT]) in the phase space The system (jM|) - (jg?)) 

admits the following critical points 

1. The critical point Pi with coordinates ip = 0,(72 = —1,(74 = 0,(75 =0 always exists. 
The eigenvalues of the linearisation around the critical point are 0, |, 1 — m (^_ 3 t) — 
7+2. This means that the critical point is non-hyperbolic thus the Hartman-Grobmann 
theorem does not apply. However, by applying the Invariant Manifold theorem, we 
obtain that: 

(a) Pi has a 1-dimensional center manifold tangent to the 99-axis provided N ^ 
and M 7^ — ^^S^ (otherwise the center manifold would be 2- or 3-dimensional). 

(b) Pi admits a 3-dimcnsional unstable manifold and a 1-dimensional center manifold 
for 

i. N < V6, < 7 < |, M > ~^Ep\ or 

ii. n < VE, I < 7 < 2, M < -^Ep-. 

In this case the center manifold of Pi acts as a local source for an open set of 
orbits in (l98l). 



(c) Pi admits a 2-dimensional unstable manifold, a 1-dimensional stable manifold 
and a 1-dimensional center if 

i. N > V6, < 7 < |, M > -^Er 1 ; or 



'For notational simplicity we will denote the image of f! e under / by the same symbol. 
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ii. N > y/6, | < 7 < 2, M < - 



Vej-y-2) , 

37-4 ' 



or 



N < VE, < 7 < |, M < -^Ep; 



3 7 -4 ' OT 

iv. TV < V6, I < 7 < 2, M > -^t 2 ^ 1 - 

(d) Pi admits a 1-dimensional unstable manifold, a 2-dimensional stable manifold 
and a 1-dimensional center manifold for 

i. N > VE, < 7 < |, m < -4^; OT 

ii. N > y/6, I < 7 < 2, M > -^EP-. 

2. The critical point P2 with coordinates 95 = 0,(72 = 1)CT4 = 0,(75 = always exists. 
The eigenvalues of the linearisation around the critical point are 0, |, 1 4- ^)=, —7 + 

M (^T 4 ) + 2. As before, let us determine conditions on the free parameters for the 
existence of center, unstable and stable manifolds for P2. 

(a) If N / — y/6 and M / there exists a 1-dimensional center manifold 
tangent to the c/?-axis, otherwise the center manifold would be 2- or 3-dimensional. 

(b) P2 has a 3-dimensional unstable manifold a a 1-dimensional center manifold (tan- 
gent the ^-axis) if 

i. N > -y/6, < 7 < |, M < or 

ii. N > -y/6, I < 7 < 2, M > ^gt_a. 

In this case the center manifold of Pi acts as a local source for an open set of 
orbits in (l98l). 



(c) Pi has a 2-dimensional unstable manifold a 1-dimensional stable and a 1-dimensional 
center manifold if 

i. N < -y/6, < 7 < f , M < or 

ii. JV < -y/6, § < 7 < 2, A/ > or 
hi. N > -y/6, < 7 < f , M > ^r 4 2) ; or 
iv. TV > ->/6, I < 7 < 2, M < 

(d) P2 has a 1-dimensional unstable manifold a 2-dimensional stable and a 1-dimensional 
center manifold if 

i. N < -y/6, < 7 < f , M > or 

ii. N < -V6, I < 7 < 2, M < ^Ep. 

3. The critical point P3 with coordinates = 0, a 2 = yg(^~2) ' Ui — 0, (J5 = exists for 

(a) 0<7<i,-^<M<^ ;O r 

(b) |<7<2,4^<A/<-feH). 
The eigenvalues of the linearisation are 

n , _ (37-4)((37-4)M 2 -2 7 +4) _ -Af 2 (4-37) 2 +6(7-2) 7 +2MAf(37-4) > _ 6( 7 -2) 2 -M 2 (4 

U ' Al — 12(7-2) ' A2 ~~ 12(7-2) ' A 3 — 12(7-2) 

As before, let us determine conditions on the free parameters for the existence of cen- 
ter, unstable and stable manifolds for P3. 
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(a) For 7, N and M such that Ai 7^ 0, A 2 7^ 0, A 3 7^ the center manifold is 1- 
dimensional and tangent to the <^-axis. Otherwise the center manifold coud be 
2-, or 3-dimensional (it is never 4-dimcnsional) . 

(b) P3 admits a 1-dimensional center manifold and a 3-dimensional stable manifold 
for 

i 0<<7<^ V^F2 M q v M 2 (4-3 7 ) 2 -6(7-2)7 , 
1. U <. 7 < 3, 73"^ <- JW U ' JV ^ 2M(3 7 -4) ' or 

ii < -7 < ± < M < JV < M 2 (4-3 7 ) 2 -6(7-2) 7 

11. U <^ 7 <^ 3, U ^ 1V1 > JV ^ 2M(3 7 -4) 

(c) In the cases 

i 0<T<^ _ ^(7-2) M _ 72y^2 ^ M 2 (4-3 7 ) 2 -6( 7 -2)7 . 
1. U ^ 7 ^ 3, 3 7 _ 4 ^ JW \/3t^4" > ^ 2M(3 7 -4) ' U1 

ii < T < — <M< ^ (7 ~ 2) N > M 2 (4-3 7 ) 2 -6(7-2) 7 . 

iii 4 <r <r 9 <r M <r n AT -> m 2 (4-3 7 ) 2 -6( 7 -2)7 . nr 

111. 3 < 7 < 2, 3 7 _ 4 < M < U, iV > 2M(3 7 -4) ' ° r 

iv. I < 7 < 2, M = 0, N e R; or 

v. §< 7 <2,0<M< - ^ 7 7 J 4 2) , N < M2(4 2M(37- 6 ir 2)7 ' the unstable mani- 
fold is 2-dimensional (hence the stable manifold and the center manifold are 
both 1-dimensional). 

vi. Otherwise, P3 has a 1-dimensional unstable manifold. Thus, it is never a 
local source since its unstable manifold is of dimension less than 3. 

The critical point R\ with coordinates (p — 0, 2 = 0, 04 = 0, 05 = 1 always exists. 
The eigenvalues of the linearisation are 

0, |, — |, I — 7. The center manifol is 1-dimcnsional and tangent to the tp-axis. The 
unstable (stable) manifold is 1-dimensional (2-dimensional) if 7 > | otherwise it is 
2-dimensional (1-dimensional). 



5. The critical point i?2 with coordinates 02 = -^-,04 = 0,05 = ^3j-4 — exists for 
< 7 < I , M 2 > ■ The eigenvalues of the linearisation are 

0> _ M + V3^(y)- 8(2 -2) ) Hf + 2 ) • Lct us dcte ~ con " 

ditions on the free parameters for the existence of center, unstable and stable manifolds 
for R 2 . 

(a) i?2 has a 3-dimensional stable manifold and a 1-dimcnsional center manifold if 



i. < 7 < f , -2^1^^ <M< -yfiyJ^N > -2M; or 
^ < 7 < lV2j^ <M< 2 y /§ y J5£,N< -2M; or 

iii. § < 7 < f , M < -V^y^, N > -2M; or 

iv. § < 7 < |,M > y/2.f^,N< -2M; or 



v. < 7 < f ,M < -2jy^,N> -2M; or 



vi. 0< 7 < |,M>2 v /| v / i ^,iV< -2M. 

(b) By reversing the sign of the last inequality, i.e., the inequality solved for N, in 
the previous six cases we obtain conditions for i? 2 having a 2-dimensional stable 
manifold, a 1-dimensional unstable manifold and a 1-dimensional center manifold. 
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6. The critical point P4 with coordinates ip = 0, 02 = — 04 = \Jl — 05 = exists 
whenever N 2 < 6. The eigenvalues of the linearisation are 

0, \ (N 2 - 6) , i (iV 2 - 4) , \N(2M + N)- \{MN + 2)7. The conditions for the exis- 
tence of stable, unstable and center manifolds is as follows. 

(a) The center manifold is 1-dimensional and the stable manifold is 3-dimensional 
provided 

i. N = 0, M £ M,7 ^ §; or 

ii. < 7 < f , -2 < N < 0, M > ; or 
hi. 0< 7 <|,0<JV<2,M< ; or 



4 



2(jV 2 -37) 



3<7<2,-2<JV<0,M< ^ ( 3 7 _ 4 / ; or 
|< 7 <2,0<iV<2,M>^5). 



(b) The stable manifold is 2-dimensional, the unstable manifold is 1-dimensional and 
the center manifold is 1-dimcnsional provided 



i. 0< 7 <4, -\/6 < N < —2, M > 



(jv 2 -3 7 ). 



or 



3 , v"^-" ^ *,±v± ^ AT(3 7 _ 4 ) i 

ii. 0< 7 <|,-2<jV<0,M< l^"-^ ' or 

hi. 0< 7 <|,0<jV<2,M> ; or 

iv. < 7 < §, 2 < N < A M < 2 ff ^ ; or 

v. § < 7 < 2, - V6 < N < -2, M < ; or 

vi. § < 7 < 2, -2 < TV < 0, M > ; or 

vii. f< 7 <2,0<iV<2,M< ; or 

viii. § < 7 < 2,2 < TV < VQ,M > 2{ N N { 1;^ ■ 

(c) The stable manifold is 1-dimensional, the unstable manifold is 2-dimensional and 
the center manifold is 1-dimcnsional provided 

i. < 7 < f , ~ VQ < N < -2, M < %^ ; or 



ii. < 7 < |,2 < N < V6,M > 



2(W 2 -37) _ 



AT(37-4) 



or 
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111. I < 7 < 2, - V6 < N < —2, M > ^ (37 _ 4 y ; or 
iv. § < 7 < 2, 2 < TV < V6, M < ^3 ■ 

7. The critical point i? 3 with coordinates tp = 0,a 2 = — ^^-,0-4 = ^p^j 1 °"5 = ^^p 4 
exists for N 2 > 4. The eigenvalues of the linearisation are 

0, 1 (- VME^MM _ 1) , 1 ( vW-i 5j v* _ ^ _ (2M+ao(3 7 -4) . Thc cond i tions for the 
existence of stable, unstable and center manifolds are as follows. 

(a) The stable manifold is 3-dimensional and the center manifold is 1-dimensional 
provided 



2(Af 2 -3 7 ) _ 
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i. 0< 7 <|,7V<-^ H ,M>-f;or 

ii. < 7 < |, -^g < N < —2, M > -f ; or 

iii. 0<7<|,2<JV< < -f ; or 

iv. 0< 7 <|,iV>^= I M<-f;or 

v. |< 7 <2,7V<-^ f ,Af<-f;or 

vi. | < 7 < 2, < N < —2, M < -f ; or 

iV. 



Vll. 



|<7<2,2<7V<^ H ,M>-f;or 
viii. § < 7 < 2, TV > M > -f . 

(b) By reversing the sign of the last inequality, i.e., the inequality solved for M, 
in the previous eight cases we obtain conditions for P3 having a 2-dimensional 
stable manifold, a 1-dimensional unstable manifold and a 1-dimensional center 
manifold. 

The critical point P5 with coordinates 



V = °' CT 2 = Af(3 7 -4)-2Af ff 4 - X 2JV+M(4-3 7 ) ' ^ = 6X1StS 

2(2M+iV) > 3M 7 , M(3 7 -4)(M(37-4)-2AT) > 6( 7 -2) 7 , and m^^^+El < 
0. The eigenvalues of the linearisation are 



12Af+6Af-3(3Af+jV)7+y / 3 % //(7 : M,jV) 3iV(7-2)+3Af (37-4) + x/3 % //(7,A-/ ,iV) (2M+iV)(3 7 -4) , 
U > 6(M(3 7 -4)-2iV) ' 6(2JV+M (4-3 7 )) ' 6JV+3M (4-37) ' Wnere 

/(7, M, iV) = 2M 3 iV(37-4) 3 +2AfA^ (AN 2 - 6 7 2 + 37-6) (3 7 -4)-M 2 (8iV 2 - 12 7 - 3) (4- 
3 7 ) 2 + 3(7 - 2) (7V 2 (9 7 - 2) - 24 7 2 ) . The stability conditions of P 5 are very compli- 
cated to display them here. Thus we must rely on numerical experimentation. We 
can obtain, however, some analytic results. For instance, there exists at least a 1- 
dimensional center manifold. The unstable manifold is always of dimension lower 
than 3. Thus the critical point is never a local source. If all the eigenvalues, apart 
form the zero one, have negative reals parts, then the center manifold of P5 acts as 
a local sink. This means that the orbits in the stable manifold approach the center 
manifold of P5 when the time goes forward. 

The critical point Pg with coordinates 



n Ve~f ^/AP(4-37) 2 + AfAT(8-67)-6(7-2)7 . , 

f = °< a2 = M(3 7 -i)-2Ar ' ff 4 = 2j V+M(4-37) ' ^ = eX1StS fof 

M(37-4)(M(3 7 -4)-2AT) > 6( 7 -2) 7 , 2(2M+iV) < 3M 7 , and 3(M ^%~ ( l^ +N) < 
0. The eigenvalues of the linearisation are the same displayed in the previous point. 
However the stability conditions are rather different (since the existence conditions are 
different from those of P5). As before, the stability conditions are very complicated to 
display them here, but similar conclusions concerning the center and unstable mani- 
fold, as for P5, are obtained. For get further information about its stability we must 
to resort to numerical experimentation. 



Physical description of the solutions and connection with observables. Let us 

now present the formalism of obtaining the physical description of a critical point, and also 
connect with the basic observables relevant for a physical discussion. These will allow us to 
describe the cosmological behavior of each critical point, in the next section. 

Firstly, around a critical point we obtain first-order expansions for P, a, <p, and p and 
p r in terms of t, considering equations: (|26|) ; the definition of the scale factor a in terms of 
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the Hubble factor H; the definition of 172; the matter conservation equations (f2~T|) and 
respectively, given by 



2H(t) = H(t) 2 (3(7 - 2)af + 3 7 (a*f + af - lj - 4a* 5 

a(t) = a(t)H(t) 
4>{t) = y/Qa*H(t) 

p(t) = —Hit) 3 (V6M(3 7 - 4K - 67) (of + of + of - l) 

p r (t) = -\2afH{tf 



(99) 



where the star-upperscript denotes the evaluation at a specific critical point. The equation 



4>{t) = -H(t) 2 ( M(3 7 - 4) ( af + of + af - 1 ) - 2 ( N a 



f + VEat 



(100) 



derived from the equation of motion for the scalar field (|29| should be used as a consistency 
test for the above procedure. Solving the differential equations and substituting the 
resulting expressions in the equation (llOOp results in 



-6M(3 7 - 4) (af + of + of - l) + 12Nof+ 
-2VQa* (37 (of + af + af - l \ - 6af - Aaf + 6^=0. 



(101) 



This integrability condition should be (at least asymptotically) fulfilled. 

Instead of apply this procedure to a generic fixed point here, we submit the reader to 
section [4] for some worked examples where this procedure has been applied. However we 
will discuss on some cosmological observables. 

We can calculate the deceleration parameter q defined as usual as [120] 



aa 



(102) 



Additionally, we can calculate the effective (total) equation-of-state parameter of the uni- 
verse w e j, defined conventionally as 



tdeff : 



Ptot 

i 

Ptot 



(103) 



where ptot and ptot are respectively the total isotropic pressure and the total energy density. 
Therefore, in terms of the auxiliary variables we have 



9 = -~(7 - 2)al - M + 1(4 -3 7 K 2 + 1(37-2) 
cff = (2-7)^-7^ + 1(4-37)^+7-1. 



(104) 
(105) 



First of all, for each critical point described in the last section we calculate the effective 
(total) equation-of-state parameter of the universe w D s using (|105[) . and the deceleration 
parameter q using (|104j) . The results are presented in Table 01 Furthermore, as usual, for 
an expanding universe q < corresponds to accelerating expansion and q > to decelerating 
expansion. 
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Table 4: Observable cosmological quantities, and physical behavior of the solutions, at the 
critical points of the cosmological system. We use the notations Mi (7) = ^ 2l< fj 3 ^ — - , 



M 2 ( 7 ) = 



4-37 



Cr.P. 


5 


Wcff 


Solution / description 


Pi 


2 


1 


Decelerating. 


P 2 


2 


1 


Decelerating. 


P 3 


-Af 2 (4-37) 2 +2 7 (37-8)+8 


-^+7-1 


Accelerating for 

< 7 < §, -Mi (7) < M < Mi (7) 


4(7-2) 






I (AT* - 2) 


I - 3) 


Accelerating for 
-V2 < n < \/2 



powerlaw- inflationary 



p 5 


3(M+N)f-2(2M+N) 


M(4-3 7 )-2JV( 7 -l) 


Accelerating for 

3(M+JV) 7 -2(2M+iV) „ n 
2AT+M(4-37) ^ U 

mattcr-kinetic-potential scaling 


2N+M(i-3-f) 


M(37-4)-27V 




3{M+N)i-2(2M+N) 


M(4-3 7 )-2JV(7-l) 


Accelerating for 

3(M+Af) 7 -2(2M+W) „ „ 
2W+M(4-37) ^ U 

matter-kinetic-potential scaling 


2N+M(4-3-y) 


M(37-4)-2W 


Ri 


1 


1 
3 


Decelerating. Radiation-dominated. 


R2 


1 


1 
3 


Decelerating, 
radiation- kinetic-potential 


R3 


1 


1 

3 


Decelerating, 
radiation- kinetic-potential scaling. 
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3.4.8 The flow as <j> -> -oo 



With the purpose of complementing the global analysis of the system it is necessary inves- 
tigate its behavior as cj) —> — oo. It is an easy task since the system (|41 | -(|45 p is invariant 
under the transformation of coordinates 



(0,cr 2 ) ->-(</», cr 2 ), V->U, X~>2, 

where U{<p) = V{—<j>) and S(0) = %(—</>). Hence, for a particular potential V, and a par- 
ticular coupling function the behavior of the solutions of the equations (|4"T)l -(|4"5 jl around 
(j> = — oo is equivalent (except for the sign of </>) to the behavior of the system near <fi = oo 
with potential and coupling functions U and S, respectively. 

If £/ and S are of class £^_ , the preceding analysis in E e can be applied (with and adequate 
choice of e) . 

The set of class C k functions well behaved in both +00 and —00 is denoted by £ k . 
Latin uppercase letters with subscripts +00 and —00, are used respectively to indicate the 
exponential order of £ k functions in +00 and in —00. 



4 Examples 

In this section we apply the mathematics discussed in previous sections to several worked 
examples from both analytical and numerical viewpoint. 

4.1 Numerical evidence of the result of theorem [8] 

For the particular case x(°i) = e MtTl , from equation ([90]). follows that er 2 = 02c '■= —r£^\ 
is an invariant set. Given | < 7 < 2, the existence conditions lead to Mi (7) < M < 

"f?, where M i(t ) = ^g^±* . For such values §£\ a2c = > 0. 

Thus the asymptotic phase configuration a\ — > — sgn<72 C oo, er 2 — > a 2 c is never approached 
(for an open set of orbits) as r — > 00. For the original system (i.e., taking the time reversal 
transformation) this means that this asymptotic phase configuration is never approached 
towards the past. In figure Q] we show the qualitative dynamics of the flow of (|84|) for the 
choice M = -y/2/3, and 7 = 1. In order to compactify the phase space we have introduce 
the coordinate transformation o\ — > tanhci. The mentioned asymptotic configuration is 
represented in figure Q] by P 3 ± . 

4.2 Quadratic gravity: F(R) = R + aR 2 . 

Quadratic gravity, F(R) = R + aR 2 , is equivalent to a non- minimally coupled scalar field 
with the potential 

V(<t>) = ^~ (l-e-V^A 2 ( 106 ) 

and coupling function 

X(4>) = eVT*. (107) 

Observe first that 

W x (4>) = x'W/xW ~ VV3 = lim W x (0) = (108) 



^+00 



and 

IT, 



= V"(^)/^(0) = - (a/873) /(l - ^^(fl = 0. (109) 
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Figure 1: Qualitative dynamics of the flow of (|84l) for the choice M = y/2/3 and 7 = 1. 
Observe that P±~ or P^ are the local attractors (sinks). P% is a local source, P 3 ~ is a saddle 
for the full dynamics, but it is a local source in the invariant set tanheri = — 1. The thick 
dashed line is an invariant set which is unstable. In fact all its points including Rf and 
Ri , act as saddle points. They correspond to cosmological radiation- dominated solutions. 



In other words, the coupling function (|107l) and the potential (I106|) are WBI of exponential 
orders M — y/2/3 and N = 0, respectively. 

It is easy to prove that the coupling function (|107[) and the potential (|106|) are at least 
E\, under the admissible coordinate transformation 

ip = r 1 = f {<!>)■ (no) 

Using the coordinate transformation (|110p we find 

W x ( v ) = o. (111) 



W V (<P) - 



and 



l-e\ 







ip > 
<p = 



(112) 



Lf > 

if = 



(113) 
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-1 -1 



tanh^) 



Figure 2: The graphic illustrates the result of theorem^ We set M — y/2/3 and 7 = 1. 
The point with tanhui = ±1 are the local sinks (thus for the original system (|41l) - (|45[) the 
scalar field almost always diverges towards the past). 



Table 5: Location of the critical points of the flow of defined in the invariant set 

{p e n e : ip = 0} for M = ^/2/3 and N = 0. 



Label 


(02,04,0-5) 


Existence 


Stability a 


Pi 


(-1,0,0) 


always 


unstable for I ^ < 7 < 3, or 
I f<7<| 
saddle, otherwise 


P2 


(1,0,0) 


always 


unstable 


P 3 


(1-35^.0.0) 


f 0< 7 <|, or 
1 I<7<1 


saddle 


Ri 


(0,0,1) 


always 


saddle 


Pi 


(0,1,0) 


always 


stable 



The stability is analysed for the flow restricted to the invariant set tp = 0. 
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In this example, the evolution equations for ip, a%, 04, and 175 are given by the equations 
(|M l) -(^7 |l with M = v/2/3 and N = 0, and W Xl W v , and J 7 , given respectively by (ITTTjl . 
(II 1 2[) and (|113p . The state space is defined by 

n e = {((p,<T 2 ,(7 4 ,a 5 ) G M 4 : < </? < e, of + a\ + a? D < l,cr 4 > 0, cr 5 > 0}. 

Let us analyse the local stability of the critical points of the corresponding system. In the 
above analysis we are not taking into account perturbations in the i^-axis. It is obvious, from 
the previous analysis, that the center manifold of these critical points contains the i/j-axis 
as a proper eigendirection. In table [5] are summarised the location, existence conditions and 
stability F°l of the critical points. 




Figure 3: Proyection of some orbits of (194p -(|97 j) in the invariant set ip = for the coupling 
function (|107[) and the potential (|106|) for 7=1. Observe that P\ and Pi are local sources, 
Ri, and P3 are saddles (P3 is the local attractor in the invariant set y — 0) and P4 (the de 
Sitter solution) is the local sink in the invariant set ip = 0. 



Let us discuss the stability properties of the critical points displayed in table [5] 

The critical point Pi always exists. Its unstable manifold is 3D provided < 7 < | or 
3 < 7 < I . Otherwise its unstable manifold is lower dimensional. 

The critical point Pi always exists. It has a 3D unstable manifold. Although non- 
hyperbolic our numerical experiments suggest that it is a local source. 

The critical point P3 exists for 0<7<|or|<7<| and it is neither a sink nor a 
local source. 

The critical point Ri always exists and it is neither a sink nor a local source. 

10 The stability is analysed for the flow restricted to the invariant set ip = 0. 
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The critical point P4 (corresponding to the de Sitter solution) always exists. Its stable 
manifold is 3D. Since P4 is nonhyperbolic the linear stability analysis is not conclusive. 
Thus we need to resort to numerical experimentation or alternatively we can use more 
sophisticated techniques such as normal forms expansion or center manifold theorem. Due 
its relevance, the full stability analysis of P4 is deserved to section [4. 2. II 

Let us discuss some physical properties of the cosmological solutions associated to the 
critical points displayed in table O 

• Pi, 2 represent kinetic-dominated cosmological solutions. They behave as stiff-like mat- 
ter. The associated cosmological solution satisfies H = 3t l , a = \/3t — C1C2, 4> — 

C3 ± y| In (3i — ci) , where Cj, j — 1, 2, 3 are integration constants. These solutions 
are associated with the local past attractors of the systems for an open set of values 
of the parameter 7. 

• P3 represents matter-kinetic scaling cosmological solutions such that H — t (3 7 ^gjZ.3( 2 |_ 2 
a = (i(3 7 -8)-3( 7 -2) Cl ) 1+ ^c 2 ,p = {t{ ^_™I$_ 2)ci y + c 3 , and cj> = c 4 + 
vj?(37— 4) in(t^37-8)-3(7-2)ci) wnere j _ ^ 2, 3, 4 are integration constants. 

• Ri represents a radiation-dominated cosmological solutions satisfying H = 2 t-c ' a = 

y/2t - CiC 2 , pr = (2t- ci )^ + C 3- 

• P4 represents a de Sitter solution with H = y~^, a = c\ exp 

In the figure [3] are are displayed typical orbits of (|M)) - (p?Tl) in the invariant set ip = 0. 
The critical points Pi and P 2 are local sources, R\, and P3 are saddles (P3 is the local 
attractor in the invariant set y = 0) and P4 (the de Sitter solution) is the local attractor in 
the invariant set ip = 0. However, concerning the full dynamics, it is locally asymptotically 
ustable as we prove in next section by explicit calculation of the center manifold at P4. 



3 L 



,V(<f>) = V . 



4.2.1 Stability analysis of the de Sitter solution in quadratic gravity 



In order to analyze the stability of de Sitter solution we can use center manifold theorem. 
Let us proceed as follows. First, in order to remove the trascendental function in Wv, let 
us introduce the new variable 

1 

u = 



1 — exp 



taking values in the range 



1 — exp 



A 



< u < 0. 



In this way we obtain the new system of ordinary differential equations 
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, 2lt 2 <72 2W(T2 

^ = 3 a"' 

-2 = ^-+(l-2J- 2 + g(37-4)a 2+( T2 (- — + g(4-3 7 )<7 5 + — j + 
+ ^(37 - 4)^ + i(37 - 4)a 5 2 + 1(4-37), 

-4 = — + (1-2) -2-4- — +a 4 ^-( 4-3 7 K + -j, 

4 = (l - 2) cxfag - 5704*6 + J(4 - 3 7 )a 5 3 + i(3 7 - 4)a 5 (114) 
describing the dynamics of quadratic gravity as <j) — > +00. 




Figure 4: Proyection of some orbits of (| 1 15[) m f/ie space 2/1,2/2,2/3 /or £/ie coupling function 
(|107|) and £/ie potential (I106|) /or 7 = 1. TTie graphic shows the behavior in the stable manifold 
of P4. TTie bulk of orbits in front of and at the right hand side of the figure represents a 
projection of the center(s) manifold(s). 



Proposition 6 The equilibrium point q := (u, cr 2 , 04, <r 5 ) = (0,0,1,0) of the system (|114[) 
is locally asymptotically unstable. 

In order to determine the local center manifold of (jl 14[) at q we have to transform the 
system into a form suitable for the application of the center manifold theorem (see section 
12.1.51 for a summary of the techniques involved in the proof). 

Proof. 

Case 7 7^ 1 



4G 




Figure 5: Proyection of some orbits of (|115[) in the space x,yi,ys for the coupling function 
(11071) and the potential (|106[) for 7 = 1. The graphic shows the unstable character of P4 
(trajectories starting at x < move away from the origin). 



3(7-1) 



- 1 



and 



Let be 7 / 1. The Jacobian of f| 1 14[) at q = (0,0,1,0) has eigenvalues 0, — 1, — -I 
and —7 with corresponding eigenvectors (0, 0, 0, 1) T , (0, 0, |, 0) T , ( 1, 0, 1, 

(0, 1, 0, 0) T . We shift the fixed point to the origin by setting 174 = 0-4 — 1. In order to trans- 
form the linear part of the vector field into Jordan canonical form, we define new variables 
(x, y%, y 2 , 2/3) = x, by the equations 



so that 



3x 

0-2 = x + 2/1 + 2/3 



3(7-1) 



CT 4 

0-5 



2/3, 
2/2, 



where 



/ x' \ 


/ 


y'1 




2/2 




V y's J 


V 


2,2/3) = 


= a; 3 + 



x yi + 3FF1) 





( 


X 


\ 


( f(x, 2/1, 2/2, 2/3) \ 






2/1 




, .9i(^,2/i,2/2,2/3) 






272 




.92(^,2/1,2/2,2/3) 


) 


\ 


273 


) 


V 33(^,2/1,2/2,2/3) / 



£ 2 2/3 



2^1 

3 



2a:yi 



2z;/3 
9(7-1) 



+ 



2xy 3 
3 : 



(115) 
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5—* 

R1 ' < 



0.4^ 




Figure 6: Proyection of some orbits of (|1 15|) in the space cc, 2/2,2/3 for the coupling function 
(|107p and the potential (|106[) /or 7=1. TTie graphic shows the unstable character of P± (the 
orbits depart from the origin for x < OJ. 



gi(x, yi ,y 2 ,y 3 ) = - fz 2 2/i7 + 2x 2 Vl + x 2 y 3l + ^-M- _ I2LM + _i_L_ _ l xy ^ + 

o„„2 I <y I gj/iya 17x3/i2/3 , 2xj/i 1 2„, , 2a;jy| aj/ 2 7 3 16a:j/ 2 7 2 157mj 2 7 

dxyj + Zxy iy3 j + — § h — 5 ^2 7 + -3 7-zrp + 3(7-1)* is( 7 -i) 2 + 

9(7-1)2 ~ -2/37 + 9 ( 7 _i) + — — +2/1+ 2/l2/37 + 3fr-TJ 3 6-&7 + 6-67^ "2^12/27 + 

2yi2/l _ j/i«37 3 , lOj/iysV _ 73i/i y| 7 , Wyiyj _ j/ij/37 3 1 i/ij/37 2 1 51/11/37 _ 161/11/3 1 vl 1 
3 (7-I) 2 T 3(7-l) 2 18(7-1) 2 "T" 9(7-l) 2 (7-I) 2 ^ (7-I) 2 " t " 3( 7 -l) 2 9( 7 -l) 2 ~T I8-I87 T 

2/1 l/jV I 5y 2 7 14j/l 

6 3( 7 -l) 3 + 6(7-l) 3 27(7-l) 3 ' 

92(x, 2/1,2/2, 2/3) = -5^2/27 + - 2 2/2 - ^2/12/27 + 2xyxy 2 + 2:2/22/37 - ^-i] + 3*7-1) ~ 2 -2/22/3 - 

±2/ 2 2/27 + 2/i2/2 + 2/i2/22/37-7^ 
37^rfj + 9(^rf)2 + 2/22/1 - 2/22/37 and 

53(2;, 2/1,2/2, 2/3) = -^ 2 2/ 3 7 - ^ - X2/12/37 + a;?/i2/3 - xy-fy + -2/i + -2/i7 - 377=77 + 3^1) ~ 
-2/1 +-2/37 " sg¥iy + s#-IJ -*» " ^?2/37 + 2/ 2 2/3 - 4 I +^+^% 2 7 - #S) + gg-fy " 
22/U/l + 2/12/37" 37^) + 37^1) " ^ ~ »7 + " 4 1 + ¥ + 3?f-TJ " ~ 

3 27/3 - ^3 . 3 . y^~y y^j 9 2 2j/f i 2/a , 2 

2/37 - 3(7-1) + 9( 7 -l) 2 + ^3 + 3(7-1) ~~ 18( 7 -1) 2 ~~ ^7 ~ 3(7-1) + 9( 7 -l) 2 + y 3 • 

The system ()1 15p is written in diagonal form 

x = Cx + f(x,y) 

y' = Py + S(x,y), (116) 
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where (x, y) € 



5 , C is the zero lxl matrix, P is a 3 x 3 matrix with negative 



eigenvalues and /, g vanish at and have vanishing derivatives at 0. The center manifold 
theorem [5] asserts that there exists a 1-dimensional invariant local center manifold W c (0) 
of pi6[) tangent to the center subspace (the y = space) at 0. Moreover, W c (0) can be 
represented as 

W c (0) = {(x,y) e E x K 3 : y = h(x) , |x| < 6} ; h (0) = 0, Dh (0) = 0, 

for S sufficiently small (see definition 2}. The restriction of (|116[) to the center manifold is 
(see definition fl3|) 

x f = f(x,h(x)). (117) 

According to Theorem [6l if the origin x = of (|117|) is stable (asymptotically stable) (un- 
stable) then the origin of l|116j) is also stable (asymptotically stable) (unstable). Therefore, 
we have to find the local center manifold, i.e., the problem reduces to the computation of 
h(x). 

Substituting y = h (x) in the second component of pi6[) and using the chain rule, 
y' = Dh (x) x', one can show that the function h (x) that defines the local center manifold 
satisfies 

Dh (x) [f (x, h (x))] - Ph (x) - S (x,h (x)) = 0. (118) 

According to Theorem [71 equation (|118[) can be solved approximately by using an approx- 
imation of h (x) by a Taylor series at x = 0. Since h (0) = and Dh (0) = 0, it is obvious 
that h (x) commences with quadratic terms. We substitute 

a 2 x 3 + O (a; 4 ) 
hx 2 + b 2 x 3 + O (x 4 ) 





hi (x) 




a\x 2 


Mz) =: 


h 2 (x) 




b\x 2 




_ h 3 (x) 




o 

C\X 



c 2 a: u 



O 



into (|118|) and set the coefficients of like powers of x equal to zero to find the unknowns 
ai,6i,ci,.... 

Since y 2 absent from the first of (|116[) . we give_only the result for hi (x) and /13 (x) . We 
find a 1 



7 



6(7-1) 



, a 2 



_ 37-5 

9(7-1) 



■J,P2 



|.0 Therefore, JTHJ yields 



2^ 
3 



5x 3 
"9" 



4x 4 
"9" 



(119) 



It is obvious that the origin x = of (|119[) is locally asymptotically unstable (saddle 
point). Hence, the origin x — of the full four-dimensional system is unstable. 
Case 7 = 1 

Let be 7 = 1. The Jacobian of (|114p at q = (0,0, 1,0) has eigenvalues — 1,-1,— |, and 

with corresponding eigenvectors (0, 0, 0, 0, |) T , (1, 0, 0, if , (0, -3, 0, 0) T and (0, 0, 1, 0) T . 
As before we shift the fixed point to the origin by setting <f 4 = (74 — 1 and define new 
variables (x, y±, y 2 , j/3) = x, by the equations 

3x 

u = — , 
2 ' 

<t 2 = x + y 1: 
04 = -3y 3 , 
o"5 = 2/2 



so that 



n We find 61 = 62 = 0. 
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/ x> \ 



















( x \ 




( /(^ 2/1,2/2,2/3) \ 


y'l 







-1 





1 






Vi 


+ 


gi(x,yi,y 2 ,y 3 ) 


y' 2 










2 

3 









V2 




92(x,yi,y 2 ,y3) 


\ y' 3 J 




V o 








-1 


J 




\ J 




\ 93(x,yi,y2,y3) j 



(120) 



where f(x, yi, y 2 , 2/3) = x 3 + x 2 yi 



2x A 

3 



2xyi 
3 



5i(z, 2/1,2/2,2/3 



2 



o 2 



■''.Vl 
3 



■'•.'/■" 
6 



9xj4 
2 



3xy 3 



£L 

6 



yry 2 

6 



„> ■ zym - f - %S 



g 2 {x,yi,y 2 ,y 3 ) = ^ + a: 2/i2/2 



S3 (a, 2/1,2/2, 2/3) = 



Vl V2 

2 



VlV3 
2 



+ S^ya and 



Ei. 

6 



V'2 f 3 

6 



18 



2 



M 

2 



Observe that the system (I120p is now in the canonical form (|116l) . Then, we proceed to the 
caculation of the center manifold. The procedure is fairly systematic and since we present 
it completely in the previous analysis we consider do not repeat it here. Instead, we present 
the relevant calculations. We obtain a\ = §,02 = = 0, b 2 = 0, c\ = \,c 2 
Taylor expansion coefficients of 



I for the 





hi (x) 




a\x 2 


h(x) =: 


h 2 (x) 




b lX 2 




_ h 3 (x) 




9 



a 2 x" 



c 2 x 



O 



By substituting this values of the unknows a%, 6i, c\, ... we obtain that, for 7=1, the dynam- 
ics of the center manifold in given also by equation (|119p . The conclusion is straighforward: 
the origin x — of (|119l) is locally asymptotically unstable (saddle point). Hence, the origin 
x = of the full four-dimensional system is unstable. 
This completes the proof. ■ 

The result of proposition [6] complements the result of the proposition discussed in |115] 
p. 5, where it was proved the local asymptotic instability of the de Sitter universe for 
positively curved FRW models with a perfect fluid matter source and a scalar field which 
arises in the conformal frame of the R + aR 2 theory. 



4.3 i? n -gravity. 

Let us consider the model with F(R) = R n where we have reescaled the usual multiplicative 
constant. It can be proved that, for n > 1, i? n -gravity is conformally equivalent to a 
non-minimally coupled scalar field with a positive potential 



V{<t>) 



= r(n)e A < n >* 

vT("-2) 



(121) 



where r(n) = \{n— l)n and A(n) = — ^li ~ , with coupling function given by (|107l) . 
In this example, the evolution equations for ip, <t 2 , <7 4 , and 175 are given by the equations 

(EU-dHJ) with M = y/2/3, N = - ^| ( "~ 2) , W x ((p) = W v {<f) = 0, and 7, given by (ITT31) . 
The state space is defined by 

£ = {{ip,<T 2 ,<7 47 a 5 ) £ R 4 : < ip < e,cr 2 



o% < 1,CT4 > 0,CT 5 > 0}. 



In table [5] and [7] are summarized the location, existence conditions and stability of the 
critical points. 
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Table 6: Location of the critical points of the flow of (|94 |) - (|97|) defined in the invariant 

V?("-2) 



set {p G tt e : if — 0} for M -- 

1(4 + y/8),N + = A (11 + 2VW) ,I» = -fi^t, and L + ( 7 ) = 



'2/3 and N - n , 

-9n+6n 2 ' 



. We use the notations 

7 2 +487 
127+4 



Label 



(02, 04, 1T5) 



Existence 



(-1,0,0) 



always 



P 2 



(1,0,0) 



always 



4 ~ 3 7 n n 

3(2-7) ' U ' U 



0<7< |, 

I <7< ' 



or 



Ri 



(0,0,1) 



always 



P 



n-2 
3(n-l) ' 



9(n-l) 2 ' ' 



n > 



2(n-l) -y/2(n-l) ■ v /(8-5Ti)n-2 
n-2 ' n-2 ' n-2 



1 < n < tia 



ft 



3(n-l) 7 
3-yn— 2n— 37 ' 



V / 2y / n^lV4n-37 
3771— 2n— 37 



,0 



0< 7 <|, Kn<r+( 7 ), or 



I <7< I 



^<«<r + ( 7 ). 



P, 



3(n-l)7 v / 2 v ^ri^T v /4n-37 



37n— 2n— 37 ' 37n— 2n— 37 



3 ^ ' — 3' '* 4 



a In this case Pg and P3 coincides. Thus, the critical points are nonhyperbolic. There exists 
a 1-dimensional stable manifold and a 1-dimensional unstable manifold provided \ < 7 < -I . 



Let us discuss the stability properties of the critical points displayed in table [7] 
The critical point Pi always exists. Its unstable manifold is 3D provided < 7 < |, n > 
I or I < 7 < 2, n > |. Otherwise its unstable manifold is lower-dimensional. 
The critical point P3 always exists and it is a local source. 

The critical point P3 exists for 0<7<|or|<7<2. It is neither a source nor a sink. 
The critical point R\ always exists and it is neither a source nor a sink. 
The critical point P4 exists for n > | . Its stable manifold is 3D if 7 ^ | , n > 2 or 
n + < ?i < 2,L(rt) < 7 < | or n + < n < 2, |< 7 <2. Where we have defined the 

grouping constants n + = -|(4 + v / 6),T(n) = — ^g^^ ■ Otherwise its stable manifold is 
lower-dimensional. 

The critical point P3 exists for 1 < 11 < n+. Thus P4 and P3 are in the same phase 
portrait for | < n < n+. P3 admits a 3D stable manifold for | < 7 < 2,N + < n < n + or 
|<7<2,l<n< N+, where we have defined N + — ^ (ll + 2%/TO) ■ Otherwise its stable 
manifold is lower-dimensional. 

The critical point P 5 exists for n = |, 7 = | or < 7 < |, 1 < n < L + (7), or 
4 < 7 < | 7 32 < n < r+ ( 7 ) ; where we have defined r+ ( 7 ) = 97+^+487+16 +4^ n admits 

a 3D stable manifold provided || < n < 2, 0< 7 < L(n), or 1 < rt < JV+, < 7 < |, or 

iV + < n < §, < 7 < ^IfeS^F " i^Xt,Tly Otherwise its stable 
manifold is lower-dimensional. 
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Table 7: Stability of the critical points of the flow of (f94 |) - (|97|) defined in the invariant 

v/j("-g) 



set {p G tt e : ip = 0} for M 



'2/3 and N 



-. We use the notations bj 



1(4 + v^),iV + = ^ (11 + 2>/I5) ,r(») = -3^,, and T + ( 7 ) = W^Wi6 + 4 



127+4 



Label 



Stability 3 



Pi 



unstable for < 7 < | , n > 



or I < 7 < |,n > f, 



saddle, otherwise 



P-2 



unstable for 7 ^ ^ , n > 1 
saddle, otherwise 



saddle 



#1 



saddle 



P 4 



stable for 

7 ^ |, n > 2 or 
h < n < 2,T(n) < 7 < I 
n+ < n < 2, § < 7 < 2. 
saddle, otherwise 



i?3 



stable for 

< 7 < 2,N + < n < n + 
|<7<2,l<n<iV + 
saddle, otherwise 



P 5 



stable for 



§ < n < 2, 0<7< r(n) or 



1 < n < N + , < 7 < 



iV+<n<|i, 0<7< 



4V96« 5 -272n 4 + 230n :i -5Ort 2 
3(4n 3 -24n 2 +29n-9) 

saddle, otherwise 



2(l0n 3 -19n 2 + 13n) 
3(4n :1 -24n 2 +29 n— 9) 



The stability is analysed for the flow restricted to the invariant set ip = 0. 



The critical point Pq exists for |<7<|,n=^.In this case Pq and P3 coincides; thus, 
the critical points are non-hyperbolic with a 2D center manifold, a ID unstable manifold 
and a ID stable manifold. 

Let us discuss some physical properties of the cosmological solutions associated to the 
critical points displayed in table [5] 

• Pi t 2 represent kinetic-dominated cosmological solutions. They behave as stiff-like mat- 
ter. The associated cosmological solution satisfies H — 3t ]_ , a — {/3t — C1C2, <\> — 

C3 ± In (3t — ci) , where Cj, j = 1, 2, 3 are integration constants. These solutions 
are associated with the local past attractors of the systems for an open set of values 
of the parameter /i. 
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R1 




Figure 7: Projection in <p = of some orbits of the flow of (|94|) - (|97[) for M = -\/2/3, N — 

— v 3 i — • We set n — 1.251. Observe that R\, P3 are in the region of physical interest. 
These are saddle points. The critical points P4 and P3 exist and are saddle points. Pi and 
Pi are local sources and P5 is a local sink. We display some orbits in the half space 04 < 
(corresponding to contracting universes) for aesthetical purposes. Observe that Pq mirrors 
the behavior of P5 . 



• P3 represents matter-dominated cosmological solutions that satisfy H = 3t7 2 2ci > a = 
(3t 7 -2 Cl )^ C2 ,p= (3t7 l 2 2ci) . +c 3 . 

• Pi represents a radiation-dominated cosmological solutions satisfying H = 2 t-c > a = 
^2t - ac 2 , p r = (2t _ 3 Cl)2 + c 3 - 

• P4 represents power-law scalar-field dominated inflationary cosmological solutions. It 
is easy to obtain the asymptotic exact solution: 

• P5 represent matter-kinetic-potential scaling solutions satisfying 

H = 3(n-2)ff:7nKr") + 37) Cl ' a = ^ n ~ 2 )T + W 2 ~ ^ + ^ <*, and 

_ _ 6(3 7 +»(-9 7 +"(67+2)-4)) , _ , y/6(n-l) log(3(n-2)t7+(n(2-3 7 H37)ci) 
P — C 3 (3(n-2)t 7 +(n(2-37)+37)ci) 2 ' ' n-2 

• P3 represent radiation- kinetic-potential scaling solutions satisfying H = 2 t-a > a = 

A /97 — „ _ „ 3(n(5n-8)+2) 1 — . V6(n-l)log(2t-ci) 

To complete the section we present in figures [JJ and [5] a numerical elaboration of the 
model under consideration. 
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°2 



Figure 8: Projection of the orbits displayed in figure^ to 05 = 0. 

4.4 Coupling functions and Potentials of exponential orders M = 
and iV = —fi 7^ 0, respectively. 



Table 8: Location of the critical points of the flow of ([94|) - ([97l) defined in the invariant set 
{pen e :tp = 0} for M — and N = -fjt. 



Label 


(<T2, 0~4, O5) 


Existence 


Stability 81 


Pi 


(-1,0,0) 


always 


unstable if /i > —y/E 


P2 


(1,0,0) 


always 


unstable if fi < -\/6 


P 3 


(0,0,0) 


always 


saddle 


Ri 


(0,0,1) 


always 


saddle 


Pi 




I? < 6 


stable fori ° ^ f ^Y* , 
saddle otherwise 


Pb,6 




/J 2 > 3 7 


stable for < 


< 7 < |, /1 2 > 37, or 
§<7<§, 3 7 <A. 2 <fg,or, 

I l<7<|,M 2 >|S 
saddle otherwise 


R 3 


( VI 2 VM 2 -4 \ 


ImI > 2 


stable if | < 7 < 2, saddle otherwise 



The stability is analysed for the flow restricted to the invariant set tp = 0. 
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As an example let us consider x> V E £ ? of exponential orders M = and N — —p, re- 
spectively. This class of potentials contains the cases investigated in |165l 1166] (there are not 
considered coupling to matter, i.e., x{4>) = 1, in the second case, for flat FRW cosmologies), 
the case investigated in [167] (for positive potentials and standard FRW dynamics), the 
example examined in [1], etc. In table [8] are summarised the location, existence conditions 
and stability of the critical points. 

Let us discuss some physical properties of the cosmological solutions associated to the 
critical points displayed in tabled] 

• -Pi, 2 represent kinetic-dominated cosmological solutions. They behave as stiff-like mat- 
ter. The associated cosmological solution satisfies H = 3t ^_ , a = {/3t — Oyc^, 4> — 

C3 ± ^/| In (3i — ci) , where Cj, j = 1, 2, 3 are integration constants. These solutions 
are associated with the local past attractors of the systems for an open set of values 
of the parameter p. 

• P 3 represents matter-dominated cosmological solutions that satisfy H 

12 

(3t 7 -2ci) 



3*7-2ei ' 

(3t 7 - 2ci) ^C 2 , p = T^r-B^ + c 3 . 



• Ri represents a radiation-dominated cosmological solutions satisfying H — 2 t- Cl > a 



y/2t - cic 2 , p r = (2 ' 3 ri p + c 3 . 



P4 represents power-law scalar-field dominated inflationary cosmological solutions. 
As t — > + the potential behaves as V ~ Vbexp[— p<f)]. Thus it is easy to obtain the 



asymptotic exact solution: H — tfl 'i_ 2ci 1 a — (fy* — 2ci) ^C2, <t> ~ ^ In 



V (tM 2 -2ci) 2 
2(6-^) 



• Ps^ represent matter-kinetic-potential scaling solutions. As before, in the limit t 
+ we obtain the asymptotic expansions: H = 3t7 ^ 2 ci ' a = (3*7 — 2ci) ~C2, </> 

1 1 r V M 2 (3t T -2ci) 2 ' 
M L 18(2-7)7 

• P 3 represent radiation-kinetic-potential scaling solutions. As before are deduced the 



following asymptotic expansions: H — — ^ — > a — V%t~~ci c 2, ~ — In 



voti 2 (2t-c 1 y 



2t-ci ' " v 

To finish this section let us re-examine the example discussed in [T] section B.l in presence 
of radiation. 



4.4.1 The dynamics near <f> = +00 for the model with Power-law coupling and 
the Albrecht-Skordis potential 

Let us consider the coupling function 

j_ 

X(<t>)= (y) Xo(0-^o)°, a >0, const., O > 0. (122) 



Observe that 

dhxx(<f>) 



^0 



dx a{(j) - (f> ) 

for all finite value of cf>. Since lnx(<^) has not stationary points thus the early time dynamics 
is associated to the limit where the scalar field diverges. 



12 The stability is analysed for the flow restricted to the invariant set ip = 0, i.e., we are not taking into 
account perturbations in the ip-axis. 
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This choice produces a coupling BD parameter given by 



3 VXo 

This types of power law couplings were investigated in [168 from the astrophysical view- 
point. For STTs without potential, the cosmological solutions for the matter domination 
era (in a Robertson- Walker metric) are a (t) oc (lnt)^- 1 )/ 3 "*!, 4>{t) oc (lnt)a. The values 
of the parameter a in concordance with the predictions of A H are a = 1, 0.33, 3 (see table 
4.2 in [168]). 

Let us consider also the Albrecht-Skordis potential given by |169j 

V(cj>) =e~i*(A +{<j>- Bf), (123) 

This particularly attractive model of quintessence is driven by a potential which intro- 
duces a small minimum to the exponential potential. Unlike previous quintessence models, 
late-time acceleration is achieved without fine tuning of the initial conditions. The authors 
argue that such potentials arise naturally in the low-energy limit of M-theory. The constant 
parameters, A and B, in the potential take values of order 1 in Planck units, so there is also 
no fine tuning of the potential (we suppose also that p=£0). They show that, regardless of 
the initial conditions, p^ scales, with p oc p^ oc t~ 2 during the radiation and matter eras, 
but leads to permanent vacuum domination and accelerated expansion after a time which 
can be close to the present. 

The extremes of of the potential (|123[) are located at ± = 1+ * Bm ^Z 1 A ^ . They are real 



if 1 > p A. The local minimum (respectively, local maximum) is located at (f> (respectively 



±V"(cf> ± ) = -2VoVl - Ap^e-V^^ 1 -^ 2 ) < 0. 

As we investigated in section [3. 4.41 the late time dynamics of the flow of (| 5 1 1) - (f54|) is 
associated with the extremes of the potential (the critical point P2 = (0,0,0)). When we 
restrict ourselves to this invariant set, we find that the critical point associated to <fi + is 
always a saddle point of the corresponding phase portrait. The critical point associated to 
<\r could be either a stable node or a stable spiral if 

8(3 + 2^ ^ 
(3 + 4/i 2 ) 2 " A 1 

or 

A< 8(3 + 2^ 
(3 + 4 M 2 ) 2 ' 

The early time dynamics of the flow of (I5T1) - (|54"|) corresponds to the limit <j> = +00. 
Observe first that 



and 



W x {4>) = dixM/xtt) = -n— TT =► I™ W x &) = (124) 

a[<p — </>o) 

wv(<p) = tyvwvfo) + 11 = ^^fejL ^+00 Wv{4>) = °- (125) 

In other words, the coupling function (|122l) and the potential (|123|) are WBI of exponential 
orders M = and N = —p, respectively. 
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It is easy to prove that Power-law coupling and the Albrecht-Skordis potential are at 
least under the admissible coordinate transformation F^l 

<p = <r* = f {<!>)■ (126) 

Using the above coordinate transformation we find 

W^M > (127) 
, (p = 

_ r 2y 2 (By 2 -l) p 

[ , <p = o 

and 

n*)={-K ; III i 129 ) 

In this example, the evolution equations for </?, 02, (T4, and (T5 are given by the equations 
(|M |) -(|9? f with M — 0, N — -(jl, and W x {<p), W v (f) = 0, and J 7 , given by (jT27| . (p8]l and 
(|129[) respectively The state space is defined by 

^ = {(^,02, 0-4, <t 5 ) € M 4 : < </? < y/e,<T% + a\ + a\ < l,er 4 > 0,cr 5 > 0}. 

In figure |9] we show some orbits in the invariant set of + °"1 + °5 ^1 f° r V 5 = for the 
model with coupling function (|122j) potential (|123j) . We select the values of the parameters: 
e = 1.00, fi = 2.00, A = 0.50, a = 0.33, B = 0.5, and <f> Q = 0. In this case P 5 is the local sink 
in this invariant set. 



5 Conclusion 

In this chapter we have extended several results about fiat FRW models in the conformal 
(Einstein) frame in scalar-tensor gravity theories including f(R) theories through conformal 
transformation) . Particularly we have considered a cosmological model based on the action 

Sef = d^x^M i(V0) 2 - l/W+xtr'WrfftVftxir 1 ^)} , (130) 

where R is the curvature scalar, cf> is the a scalar field, and V is the potential of the scalar 
field and \ is t ne coupling function. We have consider both ordinary matter described by a 
perfect fluid with equation of state p = (7 — l)p (coupled to the scalar field) and radiation 
p r in order to describe the dynamics in a cosmological epoch where matter and radiation 
coexisted. 

We have considered scalar fields with arbitrary (positive) potentials and arbitrary cou- 
pling functions from the beginning. Then, we have straightforwardly introduced mild as- 
sumptions under such functions (differentiable class, number of critical points, asymptotes, 
etc.) in order to clarify the structure of the phase space of the dynamical system. We have 
obtained several analytical results. Also, we have presented several numerical evidences that 
confirm some of these results. 

Our main results are the following. 

13 We fix here an error in formulas B6-B9 in pp. With the choice ip = <f>~ 1 the resulting barred functions 
given by B7-B9 there, are not of the desired differentiable class. 
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Figure 9: Some orbits in the invariant set erf + o\ + o~\ < 1 for the choice of (p — /or t/ie 
model with coupling function 1122(1 potential H123\) . We select the values of the parameters: 
e = 1.00, /i = 2.00, A = 0.50, a = 0.33, B = 0.5, and O = 0. 



1. We have formulated and proved proposition [3] generalizing analogous result in [119] . 
It states that if the potential V(4>) is such that the (possibly empty set) where it is 
negative is bounded and the (possibly empty) set of critical points of V(4>) is finite, 
then, the critical point 

where 4>* is & strict local minimum for V(<f>), is an asymptotically stable equilibrium 
point for the flow. From the physical viewpoint this result is relevant since it provides 
conditions for the asymptotic stability of the de Sitter solution. 

2. After the introduction of modified normalized variables, we have proved that the phase 
space has the structure of a manifold with boundaries (see propositions 21 [5]) . We have 
devised several monotonic functions for the flow of the dynamical system which allow 
for the identification of some invariant sets. 

3. We have provided normal forms for the vector field around the inflection points and 
the strict degenerate local minimum of the potential. The normal forms have been 
used to derive rates at which the critical point under investigation is approached. The 
results obtained are in agreement with the result in proposition [3] 
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4. We have proved theorem[5]which is a generalization of the related result in P~l lll3lfTT^] . 
This result state that if \{(j>) and V(<f>) are positive functions of class C 3 , such that 
X has at most a finite number of stationary points and does not tend to zero in any 
compact set of K, then, given p, an interior point of the phase space manifold, the 
scalar field, tfi, is unbounded through the past orbit 0~(p). The relevance from the 
physical viewpoint of this result is twofold. First, the inclusion of radiation in the 
cosmic budget does not influence radically the early-time behavior of the scalar field. 
This result is somewhat expected since for small scale factor a, the energy densities of 
radiation and the scalar field goes respectively as p r ~ a~ 4 , and « ^- ~ a~ 6 (the 
last approximation is supported by theorems 4.1 and 4.2 in pQ). Second, the result 
of theorem [8] makes clear that in order to investigate the generic past asymptotic 
dynamics of the flow we must scan the region of the phase space where \<f>\ — > oo. 

5. For the analysis of the system as <fi — ¥ oo we have defined a suitable change of variables 
to bring a neighborhood of = oo in a bounded set. This method was first introduced 
in |113j (see also [162) V Having constructed a dynamical system well-suited to investi- 
gate the dynamics in a neighborhood of cj> = oo, i.e., near the initial singularity, we have 
proceeded to the identification of the equilibrium points and to the analysis of their sta- 
bility by examining the signs of the eigenvalues. We have charaterized the associated 
cosmological solutions. We have obtained in this regime: radiation-dominated cosmo- 
logical solutions; power-law scalar-field dominated inflationary cosmological solutions; 
mattcr-kinctic-potential scaling solutions and radiation-kinetic-potential scaling solu- 
tions. 

6. Using the mathematical apparatus developed in the first part of the chapter, we have 
investigated the important examples of higher order gravity theories F(R) = R + aR 2 
(quadratic gravity) and F(R) = R n . In the case of quadratic gravity we have proved in 
proposition [6l by an explicit computation of the center manifold, that the equilibrium 
point corresponding to de Sitter solution is locally asymptotically unstable (saddle 
point). This result complements the result of the proposition discussed in [115] p. 
5, where it was proved the local asymptotic instability of the de Sitter universe for 
positively curved FRW models with a perfect fluid matter source and a scalar field 
which arises in the conformal frame of the R + aR 2 theory. 

7. Finally, we have investigated a general class of potentials containing the cases in- 
vestigated in 165. 1166] . In order to provide a numerical evidence for our analytical 
results for this class of models, we have re-examined the toy model with power-law 
coupling and Albrecht-Skordis potential V{4>) — e~^(A + {4>— -B) 2 ) investigated in 
[I] in presence of radiation. 
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